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3UI.^I,IARY 


A  nethod  ox"  conformal  transformation  is  developed 
that  j.iaps  an  airfoil  into  a  straight  line,  the  line 
being  chosen  es  the  extended  chord  Una  of  the  airfoil. 
The  mapping  is  accomplished  by  operating  directly  with 
the  airfoii  ordinates.    The  absence  of  any  preliminary 
transformation  is  found  to  shorten  the  v/ork  substantially 
over  that  of  previous  methods.    Use  is  made  of  the 
superposition  of  solutions  to  obtain  a  rigorous  counter- 
part of  the  approximate  methods  of  thln-airfoll  theory. 
The  method  is  applied  to  the  solution  of  the  direct  and 
inverse  problems  for  arbitrary  airfoils  and  pressure 
distributions.    Numerical  examples  are  given.    Appli- 
cations to  more  general  types  of  regions,  in  particular 
to  biplanes  and  to  cascades  of  airfoils,  are  indicated. 


INTRODUCTION 


In  an  attempt  to  set  up  an  efficient  numerical  method 
for  finding  the  potential  flow  through  an  arbitrary  cas- 
cade of  airfoils  (reference  1)  a  method  of  conformal 
transformation  was  developed  that  v/as  found  to  apply  to 
advantage  in  the  case  of  isolated  airfoils. 

The  method  consists  in  transforming  the  isolated 
airfoil  directly  to  a  straight  line,  namely,  the  extended 
chord  line  of  the  airfoil.   The  absence  of  the  hitherto 
usual  preliminary  transformation  of  the  airfoil  into  a 
near  circle  m.akes  for  a  decided  sim.plif ication  of  concept 
and  procedure. 
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The  exposition  of  the  methocl,  followed  by  its  appli- 
cation to  the  direct  problem  of  the  conformal  mapping  of 
given  airfoils,  is  given  in  part  I  of  this  paper.   In 
part  IT  the  method  is  applied  to  the  inverse  problem  of 
airfoil  theory;  namely,  the  derivation  of  an  airfoil  sec- 
tion to  satisfy  a  prescribed  velocity  distribution.   A 
comparison  vi^ith  previous  inverse  methods  is  made.   Addi- 
tional material  that  vv-ill  be  of  use  In  the  application  of 
the  method  is  given  in  the  appendixes.   In  appendix  A  cer- 
tain numerical  details  of  the  calculations  are  discussed. 
In  appendix  B  extensions  of  the  rjiethod  to  the  conformal 
mapping  of  other  types  of  regions  are  indicated.   The 
relation  of  the  methods  used  for  the  mapping  of  airfoils 
to  the  Cauchy  Integral  formula  Is  dlscixssed  in  appendix  C. 

Acknowledgment  is  made  to  ?Trs.  T,ois  Evans  Doran  of 
the  com_puting  staff  of  the  Langley  full-scale  tunnel  for 
her  assistance  in  miaking  the  calculations. 


S^ffiEOLS 

z  =  X  +  iy   plane  of  airfoil 

t    =  ^  +  li'i   plane  of  straight  lines 

p  plane  of  unit  circle 

cp  central  angle  of  circle 

Ax      component  of  Cartesian  mapping  function  (CMF) 
parallel  to  chord 

Ay      component  of  Cartesian  mapping  function  perpen- 
dicular to  chord 

AXq,  Ay^   particular  CI/lF's,  tables  I  and  II 

T       displacement  constant  for  locating  airfoil 

r  =  2r  diameter  of  circle,  semllength  of  straight  line 

^n  ~  ^n  +  i'bn  coefficients  of  series  for  C!!?' 

Sjq-      negative  of  central  angle  of  circle,  corresponding 
to  leading  edge  of  airfoil 
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prp      central  angle  of  circle  minus  luO  ,  corresponding 
to  trailing  edge  of  airfoil 

c       airfoil  chord 

C7       section  lift  coefficient 

Vg      velocity  at  surface  of  airfoil,  fraction  of  free- 
streaia  velocity 

v>-,      velocity  at  surface  of  circle,  fraction  of  free- 


■P 


V  J'    C^  U   O  LA  -L  J.  O'.  v^  t;   V^  X    ^^  J-  J.  W  A.  ^  ^ 


strearr!  velocity 

V  free -stream  velocity 

ds  element  of  length  on  airfoil 

r  circulation 

u-t  thickness  factor 

Uq  camber  factor 

T  thickness  ratio 

X  ^  normalising  constant 

Ic  denom.inator  of  equation  (I7) 

C  camber,  percent 

6x,  5y   incremental  C^JP's 

U  positive  area  under  approximate   Vq(cp)   curve 

L  negative  area  under  approxim.ate   v-q(cp)   curve 

a  angle  of  attack 

aj  ideal  angle  of  attack 

•Y  =  ex  +  p^^ 

^,  true  potential 

^Q_  approximate  potential 

9  central  angle  of  near  circle 

c  =  CD  -  e 
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Subscripts: 

N       leading  edge  (nose) 

T       trailing  edge 

c        camber 

t       thickness 

o,  1,  2   saccessive  approxlraation  in  direct  or  inverse 
CIIF  nethods 

I  -  THE  DIRECT  POTEITTIAL  PRGBL-"2M  OF   AIRFOIL  THEORY 

THE  CARTESIAN  MAPPING  PUNCTION 
The  Derivation  of  the  Cartesian  Mapping  Function 

Consider  the  transformation  of  an  airfoil,  z-plane, 
into  a  straight  line,  t-plane  (fig.  1).    The  vector 
distance  between  conf orijially  corresponding  points  such 
as   Pg   and   Pf  on   the  tv/o  contours  is  composed  of  a 
horizontal  displacenent   Ax   and  a  vertical  displace- 
ment  Ay.   The  quantity  Ax  +  i  Lj      is  only  another  way 
of  writing  the  analytic  function  z    -    t;    that  is, 

z    -    t    =    {x   +    ij)    ~    {I   +   It,) 

=  (x  -  J)  +  l(y  -  n) 

=  Ax  +  i  Ay  (1) 

By  Rlemann ' s  basic  existence  theorem  on   conformal 
mapping,  the  function   z  -  t   connecting  conformallv 
corresponding  points  in  the  z-  and  k-planes  is  a  regular 
function  of  either   z   or  t      everjrvirhere  outside  the 
airfoil  or  straight  line.   This  f miction  will  be  referred 
to  as  a  Cartesian  mapping  function,  or  CMF.   In  order  to 
map  an  airfoil  onto  a  straight  line,  the  airfoil  ordi- 
nate s   Ay   are  regarded  as  the  imaginary  part  of  an 
analytic  function  on  the  straigh^t  line  and  the  problem^ 
reduces  to  the  calculation  of  the  real  part   Ax. 
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The  calcultition  of  the  real  part  of  an  analytic 
function  on  a  closed  contour  from  the  known  values  of 
the  imaginary  part  is  well  kno'wn.    It  is  convenient  for 
this  calculation  to  consider  the  straight  line  as  con- 
forraally  related  to  a  circle,  p-plane,  by  the  ffuniliar 
transformation 

I    ,    r    =   V   +     —  (2a) 

P 

v/here  the  constant  displacement   t   has  been  inserted  for 
future  convenience  in  locating  the  airfoil.   ?or  corre- 
sponding points  on  the  straight  line  and  the  circle,, 
equation  (2a)  reduces  to 

£=:T  +  r>COSO 

;2b) 

r,  =  0  : 

Considered  as  a  function  of   p,   therefore,  the  CMF  z    -    t 
is  regular  everywhere  outside  the  circle  and  is  therefore 
expressible  by  the  inverse  pov/er  series: 


z  -  £  =  ZI  —  (5) 

0  p"^ 

The  analogy  of  equation  (J)  with  the  Theodorsen-Garrick 
transformation  (reference  2) 

log  El :.  >~  rii 

Pip- 

v;hich  relates  confornally  a  near  circle,  p' -plane,  to  a 
circle,  p-plane,  may  be  noted.   On  the  circle  proper, 

icp 
where   p  =  Re   ,   and  defining   Cj^^^  e  a^^  +  ibj^,   equa- 
tion (5)  reduces  to  two  conjugate  Fourier  series  for  the 
ClilF;    namely, 

^^   =  a   +  yZ  -^   COS  nc?  +  Y^  —   sin  nP         (i;) 

1   R-""'  1   R 
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=   1^0    ^^ 


Ih. 


R 


n 


cos   no   -  / 


1      R- 


13: 

n 


sm  ncp 


(5) 


These  series  evidently  determine   Ax   fron  Lj      or  vice 
versa . 

An  alternative  method  of  performing  this  calculation 
is  possible.   It  is  knovvn  that  if  the  real  and  imaginary 
parts  of  a  function  are  given  by  conjugate  Fourier  series, 
as  in  equations  (Ii)  and  (5)^  with  the  constant  terms 
zero,  two  integral  relations  are  satisfied.   (See,  for 
example,  references  2  and  5;  also,  appendix  C.)   These 
relations  are 


Ax(cp) 


_1_ 
2ti 


2rr 


Ay(cr'  )  cot  - 


do' 


(6) 


Ay(cp) 


_  _1_ 
2v 


t 

Jo 


2tt 


'S>  ' 


cot: 


^ .-^  dcp' 


(7) 


Before  the  detailed  application  of  the  CT-LP  z    -    [> 
to  the  solution  of  the  direct  and  inverse  problems  of 
airfoil  theory  is  made,  some  necessary  basic  properties 
of  this  function  v/ill  be  discussed. 


Airfoil  Position  for  G-iven  CLIP 

It  is  noted  first  that  the  regions  at  infinity  in 
the  three  planes  are  the  same  except  for  a  trivial  and 
arbitrarv  translation;   that  is,  by  equations  (1),  (2a), 
and  (5),^ 


-  ^  E  Ax^  +  i  Ay, 


1  im  z 
z,  ^-^  ^ 

lim  t  =  p  +  T 

i,  p->-w 


L   +  ih 
■■0    -^^o 


>   (8) 


Secondly,  if  an  airfoil  is  to  be  mapped  into  a 
straight  line,  it  becomes  necessary  to  knov/  the  point  on 
the  straight  line  corresponding  to  the  trailing  edge  of 
the  airfoil.    For  a  given  CI/IF,   £ix(c),   Ay(<^),   and 
straight  line  of  length   2r   located  as  in  figure  1, 
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the  airfoil  coordinates 
tions  (1)  and  (2b)  as 


are    obtained  from  ecma- 


T    +   r   cos    c,i   +   Ax(cr) 


y  =  Ay(¥) 


(9) 
(10) 


The  leading  and  trailing  edges  of  the  airfoil  will  be 
taken  as  the  points  corresponding  to  the  extreraities  oi 
the  airfoil  abscissas.   The  corresponding  locations  on 
the  circle  are  therefore  determined  by  maxinizing  x 
with  respect  to   C"   in  equation  (9)'   Thus 


or 


dx 

dcf 


sm 


+ 


sin  Q^  = 


dAx 
r   dq^ 


d&x 
dcp 


11) 


?he  condition  (11)  yields  (usually  by  graphical  deter- 

corrssoondmg  to  the  leading  and 
1) 


mination)  the  angles 
edges  (fig". 


trailing 


=  -8 


-  TT 


N 


J 


(12) 


It  will  be  found  convenient  to  so  alter  the  position 
and  scale  of  a  derived  airfoil  tiiat,  for  examiole,  its 
chordwise  extremities  are  located  at  x  =  il   and  the 
trailing  edge  has  the  ordinate   y  =  0   (to  be  referred 
to  as  the  normal  f orn) .    The  chord   c   of  a  derived  air- 
foil is  by  definition  the  difference  in  airfoil  abscissa 

and  (9), 


extremities, 
c  =  r  ^cos  3- 


N 


or  by  equations  (12) 


cos 


Th£ 


increase    in   scale   ii^ori 


to 


obtained   simply  by  m.ultipl^ring 


Ax(^cnrj) 

some  desired   c 
i\x,   and  Ay 


(15) 


o 


,-1,.  ii^-  ,  v<J 

factor   Cq/c.   The  translation  necessary  to  bring  the 
trailing  edge  of  the  airfoil  to  its  desired  location 
then  accomplished  by  adjusting  the  translation  consta 
and  bQ . 


-S 

the 
is 

iitS  T 
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Velocity  Distribution  on  Airfoil 

Once  the  CMP  Ax(cp),  Ay(cp  )   and  the  diameter  of 
circle   r   of  an  airfoil  have  been  deterr.iined,  the 
velocity  at  a  point  on  its  surface  Is  obtained  in  a  v/ell- 
known  manner  as  the  product  of  the  known  velocity  at  the 
corresponding  point  of  the  circle  and  the  stretching 
factor  from  the  circle  to  the  airfoil;  that  is, 

vz(cr)  =  p  g  vp(c?)  (ll^) 

v;here   v-p(cp)   is  half  the  velocity  on  the  circle  (since 
r  =  2R)   and   ds   is  the  eleiient  of  length  on  the  airfoil. 

The  velocity  on  the  circle  v.^i'^)  ^    v;hich  makes  the 
point   cp  =  IT  +  Prp   corresponding  to  the  trailing  edge 
of  the  airfoil  a  stagnation  point  (Kutta  condition),  is 

Vp(cp)  =   sin  (cp  +  q)  +  sin  (a  +  3^)  |  (15) 

v/here   a   is  the  angle  of  attack.   The  velocities   Vp 
and   V2   are  expressed  nondimerisionally   as  fractions  of 
free-stream  velocitjr.    The  stretching  factor   ds/dcp   is 
obtained  from  equations  (9)  and  (10)  as 


ds  _  //dx\2   /dy \2    //di X       .   \2  ^  fdAy Y  t,r, 

The  velocity   V2(cp),  equation  (14),  therefore  becomes 

I  sin  ( cp  +  a)  +  sin  (a  +  S^^^l 
V2(cf')  =   ,  ^:=rz=:^^      '^^  (17) 

J         \v   dcp/ 

This  equation  is  the  general  expression,  in  terns  of  the 
CM?,  for  the  velocity  at  the  surface,  equations  (9)  and 
(10),  of  an  arbitrary  airfoil.   The  denominator  depends 
only  on  the  airfoil  geometry,  v/hlls  the  mwierator  depends 
also  on  the  angle  of  attack.    Equation  (1?)  is  similar 
to  the  corresponding  expression  in  the  Theodorsen-Garrlck 
method  except  for  the  absence  of  the  factor  representing 
a  preliminary  transformation  from  the  airfoil  to  a  near 
circle . 
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The  expi'essions  for  the  lift  coefficient  and  ideal 
angle  of  attack  nay  be  noted.   The  circulation  T   around 
the  aii^foil  is  (V   is  free -stream  velocity) 

r  =  i^n  RV  nin  (a  +  3,^.)  (l8) 

The  lift  coefficient  c-^      is  defined  by 

-  CC7V  E  r 

2    ^ 
Hence 

C7  =  l,iT  -  sin  fa  +    3,^^  ( lo ) 

G        \        -  / 

where  the  airfoil  chord   c   is  given  by  equation  ( 15 ) • 

The  ideal  angle  o£   attack  {.L^efevsnce   2)    is  defined 
as  that  angle  of  attack  for  v;hioA  a  stagnation  point 
exists  at  the  leading  edge;  that  is,   Vg  =  0   for   cp  =  -S, 
in  equation  (17)  •   Heiice, 

Ct  -  -^—  (20) 


Superposition  of  Solutions 

The  sun  of  two  analytic  functions  is  an  analytic 
function;  therefore,  for  a  given  p -plane  circle,  the  sun 
of  two  Cr.ys  is  itself  a  CMb'  as  is  also  evident  frcn 
eqiiations  {'l\.)    to  (7).   Thus,  the  CI.'F's   /Ix-i  +  i  A^y-   and 

AX2  +  i  Ay2   cif  tv7o  component  airfoils  may,  for  the  same  r, 
be  added  together  to  give  a  C?/!F   Cax]_  +  Ax2'\  ■*■  iCAy]_  +  ay2') 

and  thence,  by  equation  (17)>  an  exact  velocity  distribu- 
tion for  a  resultant  airfoil.   The  resultant  profile  and 
its  velocity  distribution  is  a  superposition  in  this 
sense  of  the  conp'^nent  profiles  and  velocity  distributions. 
Thus,  v/ithout  sacrifice  of  exactness  and  with  no  great 
increase  of  labor,  airfoils  nay  be  analyzed  and  synthe- 
sized in  terms  of  component  symrrietrlcal  thickness  distri- 
butions and  nes.n  canber  lines.   This  result  provides  a 
rigorous  counterpart  of  the  well-known  approximate  super- 
position m.ethods  of  thin-airfoil  v^'rtex  and  source-sink 
potential  theorv. 
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As  a  particular  cass  of  superposition,  a  known  Clip 
Ax  +  1  Ay  may  be  multiplied  by  a  constant   S   and  the 
resulting  ClvJP   S  Ax  +  IS  Ay   determines  a  new  profile  by 
the  new  displacements   S  Ax,  S  Ay  from,  points  on  the 
original  straight  line.    It  is  evident  that,  except  for 
the  corrections   (S  -  1)  Ax   to  the  airfoil  abscissas, 
this  nev\r  profile  is  Increased  in  thickness  and  camber 
over  the  original  profile  by  the  factor   S.    The  effect 
on  the  velocity  distribution  is  that  of  multiplying  the 
derivatives  in  equation  (1?)  by   S.   By  virtue  of  a  reduc- 
tion in  scale  by  the  factor   l/S   this  profile  ma7,'  also 
be  regarded  as  obtained  from  the  original  one  by  using 
the  same   Ax,  Ay  but  a  length  of  line   l/S   tines  the 
length  of  the  original  one . 

The  use  of  superposition  as  v^fell  as  the  application 
of  the  CMF  to  some  particular  airfoils  will  be  illustrated 
next . 


Application  oi    the  CI;1P  to  Some  Particular  Airfoils 

Symmetrical  thickness  distributions.-  The  Cartesian 
mapping  function  was  calculated  for  a  sjTnmetrical  JO- 
percent  thickness  ratio  Joakowskl  profile  from  the  known 
conformal  correspondence  between  a  Joukowski  profile  and 
a  straight  line.   The  CMF  is  given  in  normal  form  in 
table  I.   The  associated  constants   t^   and   r^   are 

given  in  table  II  and  the  profile  itself,  as  determined 
either  from  the  standard  formulas  or  from  equations  (9) 
and  (10),  is  shov/n  in  figure  2(a).   The  syrametrj'^  of  the 
profile  required  only  the  calculation  of   Ax('P),  Ay((?) 
for  0   '^   <V   ±   iQO^ .      The  corresponding  velocity  distri- 
bution (fig.  2(b))  was  obtained  from,  equation  (17)  by 
use  of  the  computed  values  of  the  derivatives .   At  the 
cusped  trailing  edge  the  velocity  as  given  by  equation  (17) 
is  Indeteriiiinate  •  ho'.vever,  the  limiting  form  of  equa- 
tion (17),  determined  by  differentiation  of  numerator 
and  denominator,  is 

icos  (cp  +  a)| 
lim  V  =  ' (21) 


r  dP' 
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It  Is  seen   from  this  expression  that  the  velocity  at  a 
Gusped  edge  depends  on  the  secoxid  derivatives  of  the 
napping  function,  that  is,  on  the  curvature  at  the  cusp. 

The  computed  second  derivatives   d  Ax^-(-/acp' 

of  the  CMF  of  table  I  are  plotted  in  figure  3  ^or   a 

of  values  of   cp   near  l80°. 


d^^y^t/d^^ 


The  Gi;!r''s  for  syriirietrical  profiles  of  different 
thickness  ratios  v/ere  determined  frora  that  for  the 
Joukowskl  profile  as  indicated  previously  in  the  section 


t5o 


Superposition  of  Solutions 


The    factor 


to   raultiply      Ax 
ratio      T 


O' 


":/o 


to  obtain 


.3  obtained  from 


u 


t 


by  which 


profile  of  thickness 


+  Uh 


ut 


^max 


Ax(c?,t) 


.x(c^tII 


-=-  =  T 


where 


.■-■o 


is  the  Kiaxim'iim  airfoil  ordinate  of  the  knov/n 


CMF  (table  I)  and  the  denominator  represents  the  seraichord 
of  the  derived  profile.   The  solution  for  uj-   is 


U+-  = 


^yc 


ma: 


r^T 


+  T  (  — ^— 
V     2 


-  Ax. 


(22) 


Values  oJ 


u 


t 


v;ere  calculated  from  this  formula  for 


thickness  ratios  of  2)4.  percent  and  12  percent  and  are 
given  in  table  II.   Eie  resulting  CMF's  were  then  nor- 


malized as  indicated 


the  section  "Airfoil  Position 


for  Given  CTiP"''  so  that  the  actual  factors 


multiply  the  original   Ax 


o^ 


Ay 


o 


we: 


\\i 


by  Vvhich  to 


t 


These  values 


are  given  in  table  II,  together  with  the  associated 
constants   t   and   r.   The  profiles  thus  determined  are 
shown  in  figure  2(a)  and  the  corresponding  x'elocity  dis- 
tributions in  figure  2(b). 

The  derived  profiles  are  not  Joukowski  profiles. 
The  point  of  m.aximum  thickness  is  shifted  back  along  the 
chord  somevd'iat  as  the  thickr^ess  ratio  decreases.   Con- 
versely, the  point  of  maximuir.i  thickness  virould  be  shifted 
forward  by  going  from  a  thin  Joukowskl  profile  to  a 
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thicker  one.   (This  result  was  the  reason  for  starting 
from  a  thick  section.)   The  CIvIP  for  the  12-percent  thick 
derived  profile  is  illustrated  in  figure  II.   It  is  to  be 
noted  that  the  horizontal  displacenent  function   AZq^(cp) 
is  symraetrical  abo^it   c?  -  tt,   whereas  the  vertical  dis- 
placement function  Ay^-j-^lcp)   is  antisyrmiietrical  about 

cp  =  TT. 

Mean  cainber  lines.-  The  CMP  was  next  calculated  for 
a  clrc'uTar-arc  profile  of  6-percent  camber  from  the  known 
conformal  cox^respondence  between  a  circular  arc  and  a 
straight  line.    The  norinpli^ed  Ci'&''  and  its  derivatives 
are  given  in  table  III.   The  CIIP  is  illustrated  in  fig- 
ure I|.   The  s^iHTiraetry  in  this  case  is  v;ith  respect  to 
CO  =   90°   and  '  c?  =  270^,   the   Ax,^^(cp)   being  antisynnetrica'J 

and  AyQ(3(a")   s^'rmnetricai.   The  circular-arc  mean  camber 
line  is  shown  in  figure  5(s)  3-^^^  the  corresponding 
velocity  distribution  in  figure  j\h) , 

Derived  mean  camber  lines  were  obtained  from  the  CHP 
for  trie  circular  arc  in  a  ruanner  similar  to  that  for  the 
symmetrical  profiles.   The  expression  determining  tho 
factor  u^   for  a  desired  percent  camber   C   is 


^c  ^Jo 

max 


2  [Tq  cos  cpjj  +  u^   Ax 


{"d 


=  C 


with  the  solution  for  u^ 


20  r^  cos  cp-^T 
Uc  -  2 S (23) 

Ay^     -  2C  Ax/co  ) 
^max         V  '■■/ 

The  angle   cpj^   in  equation  (25)  (as  in  equation  (22)) 

corresponds  to  the  extremity  of  the  derived  mean  line. 
Because  the  factor  u^   is  to  multiply  the  derivative 

dAxQ(  cp)ydcp,   the  angle   cp   as  determined  by  the  maxi- 
mum condition  (11)  depends  on  Up .   One  op  tv.o  trials 
are  sufficient  to  determine   u^   simultaneously  with   cp,j. 

from  equations  (23)  and  (11)  for  a  given  desired  camber  C. 

Values  of   u^   and   cp,.   (also   Cftyi   by  svmmetry)  are 
u         J,]  ^         ^       ^  .J 

given  in  table  IV  for  derived  cambers  of  5  and  9  percent. 

The  actual  multiplying  factor  to  obtain  the  derived 

CMP ' s  in  normal  form  is  given  in  table  IV  as  Xu^. 
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The  derived  camber  lines  are  shown  in  figure  5(2-)  • 
It  is  seen  that  the  derived  camber  lines  have  been 
separated  into  distinct  upper  and  lower  surfaces.   Fur- 
thermore, for  the  9~pe^cen"t  camber  line  the  "lower" 
surface,  that  is,  the  surface  corresponding  to  the  lower 
part  of  the  straight  line  or  circle,  lies  above  the 
"upper"  surface.   Although  such  a  camber  line  is  physi- 
cally meaningless  by  Itself,  nevertheless  its  CMF  can  be 
compounded  v/lth  that  for  a  thickness  distribution  to  give 
a  physically  real  result  (if  the  resultant  profile  is  a 
real  one).   The  velocity  distribution  of  the  J-pei-^cent 
camber  line  is  given  in  figure  '^(b).      The  "velocity  dis- 
tribution" of  the  9-pei^cent  camber  line  is  included  in 
figure  5(t')  for  aritliriietical  comparison  although  it  is 
physically  jneaningless  for  the  reason  just  mentioned. 

The  velocities  at  the  cusped  extremities  of  the 
camiber  lines  are  given  by  equation  (21).  The  second 
derivatives  of  the  ClIF   of  table  III  were  computed.   They 

are  plotted  in  figure  3  as  d^Ax^^/d(p    ,    d  Lj^^/dcp^      for 

a  range  of   cp  near  l80°.    These  second  derivatives,  in 
combination  with  those  for  the  syiraiietrlcal  profile,  can 
be  used  to  give  a  more  accurate  determination  of  the 
velocity  at  and  near  a  cusped  trailing  edge  than  is 
obtained  by  using  equation  (17)  near  the  trailing  edge. 

Combination  of  sATg-.ietrical  profile  and  mean  camber 
line . -  The  CMF's  der ived  for  the  symmetrical  profiles  and 
for  the  mean  camber  lines  can  now  be  Gom_bined  in  varying 
proportions  to  produce  airfoils  having  both  thickness 
and  cariiber.   These  airfoils  may  be  useful  in  themselves 
or,  as  in  the  following  sections,  may  be  used  as  initial 
approximations  in  both  th^e  di.rect  and  inverse  processes. 

As  an  illustration  of  such  combinations,  the  CMP 
of  the  12-percent  thick  symjiietrlcal  profile  of  fi^^ure  2(a) 
and  the  C?/IF  of  the  6-percent  camber  circular  arc  of 
figure  5(3-)  were  added  together.   The  airfoil  profile 
thus  determined  is  shown  in  figure  6(a).   For  comparison, 
the  airfoil  obtained  in  the  manner  of  thin-airfoil  theory 
(see,  for  example,  reference  ix)  by  superposition  of  the 
sajiie  symmetrical  profile  and  a  b. 5 -percent  camber  cir- 
cular arc  (in  order  to  duplicate  the  camber  of  the  exact 
airfoil  more  closely)  is  indicated  in  the  figure.   The 
velocity  distribution  of  the  dotted  airfoil  should, 
according  to  thin-airfoil  theory,  be  the  sum  of  the 
svmm.etrical-prof ile  velocitv  and  the  increment  above  the 
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free-stream  value  of  the  camber-line  velocity.   This 
velocity  distribution,  determined  from  the  tv;o  com.ponent 
exact  distributions  at  zero  angle  of  attack,  is  shovi/n 
dotted  in  figure  6(b).   The  exact  velocity  distribution 
of  the  "exact"  airfoil  of  figure  6(a)  was  determined 
for  the  same  lift  coefficient   (c-^  =  0.88,  a  =  1°15 '  ) 

from  the  known  Clip.   This  distribution  is  shown  in 
figure  6(b).   The  tv;o  velocity  distributions  differ  ap- 
preciably, although  in  the  directions  to  be  expected 
from  the  differences  in  shape  of  the  corresponding  air- 
foils. 

It  appears  that  the  CI'iF's  of  a  relatively  small 
nimiber  of  useful  thickness  distributions  and  camber  lines 
would  suffice  to  yield  a  large  number  of  useful  combi- 
nations of  vi^hich  the  (perfect  fluid)  characteristics  could 
be  determined  exactly  and  easily  in  the  manner  indicated. 

The  superposition  of  solutions  can  also  be  used  with 
the  airfoil  mapping  methods  based  on  the  conformal  trans- 
formation of  a  near  circle  to  a  circle.    There  is  a 
decided  advantage,  hov/ever,  in  working  with  the  airfoil 
ordinates  directly,  both  in  the  facility  of  the  calcula- 
tions and  in  the  insight  that  is  maintained  of  the  rela- 
tionship betv.-een  an  airfoil  and  its  velocity  distribution. 


THE  DIRECT  POTENTIAL  PROBLEM  FOR  AIRFOILS 


The  direct  problem,  for  airfoils  is  that  of  finding 
the  potential  flow  past  a  given  arbitrary  airfoil  section 
situated  in  a  uniform  free  stream.   This  problem  can  be 
solved  by  a  CMP  method  of  successive  approximation  some- 
v/hat  similar  to  that  in  reference  2. 


Method  of  Solution 

Suppose  an  airfoil  to  be  given  as  in  figure  6(a). 
The  chord  is  taken  as  any  straight  line  such  that  perpen- 
diculars drawn  from  Its  extremities  are  tangent  to  the 
airfoil.   For  examxple,  the  "longest-line"  chord,  that  is, 
the  longest  line  that  can  be  drav/n  within  the  airfoil, 
satisfies  this  definition.   The  x-axis  is  taken  along 
this  chord  and  the  origin  is  talcen  at  its  midpoint. 
Suppose,  in  addition,  an  initial  CMF  LXq      and  i^Jo, 
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straight  line  Vq,      and  chordv/ise  translation  constant   t^ 
to  be  given  such  that  the  corresponding  airfoil  has  the 
same  chord  and  is  similar  in  shape  to  the  given  airfoil. 
(At  the  worst  the  initial  airfoil  could  be  the  given 
chord  line  itself.) 

At  the  chordwise  locations   Xq(cp)   of  the  initial 
airfoil,  corresponding  to  an  evenly  spaced  set  of   up- 
values by  equation  (9),  the  differences   5y-|('P)   betv/een 

the  ordinates   Ay-j_(<P)   of  the  given  airfoil  and   Ay  ( cp) 

of  the  initial  airfoil  are  measured.   The  ordinate  dif- 
ferences  5yn(cp)   determine  a  conjugate  set  of  abscissa 

corrections   5x-j_(cp)   in  accordance  either  v/ith  equa- 
tions (ii)  and  (5)  or  equation  (6).   The  details  of  this 
calculation  are  gi\'en  in  appendix  A. 

The  initial  semilength  of  straight  line   r^   corre- 
sponding to  the  initial  airfoil  is  tlien  corrected  to   r-,  , 
and  the  translation  constant  T^   adjusted  to   t^^   go 
that  the  use  of   r-|_   with  the  first  approximate  C!-IP 
Ax-|  =   AXq  +  5x-]_,  Ljj_   ~  IxJq   +   5y-i   yields  a  first  approxi- 
mate airfoil  of  vCnich  the  chordwise  extremities  coincide 
with  those  of  the  given  airfoil.   This  correction  is 
described  in  detail  presently.   If  the  first  approximate 
airfoil  is  not  satisfactorily  close  to  the  given  airfoil, 
the  procedure  is  repeated  for  a  second  approximate  air- 
foil, and  so  on.    The  successive  airfoils  thus  deter- 
mined provide  a  very  useful  criterion  of  convergence  to 
the  final  solution;  namely,  the  given  airfoil.   Evidently, 
the  fundamental  relation  betv/een  an  airfoil  and  its 
mapping  circle 

°1    ^2 
z-p  =  Co+  —  +  -—  +.  .  . 

can  be  used  in  the  manner  indicated  to  effect  directly 
the  transformation  of  an  airfoil  into  a  circle.   It 
appears  preferable,  however,  to  subtract   R^/p   from  the 
second  term  on  the  right  and  thence  to  Introduce  the 

R 

straight-line  variable   t  =^  P  +  -5". 

The  exact  velocity  distribution  of  any  of  the 
"approximate"  airfoils  (hence  the  approximate  velocity 
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distribution  of  the  given  airfoil)  may  be  obtained  from 
equation  (I?)  using  the  derivatives  of  the  corresponding 
CMF.   The  zero-lift  angle   Prp   to  be  used  in  equation  (17) 
is  determined  for  each  approximate  airfoil  along  with  the 
corresponding  correction.  foi>   r. 

The  correction  for   r   is  necessary  because  if  the 
chordwise  locations  of  the  first  approximate  airfoil  were  ■ 
computed  by  equation  (9)  v/ith  tiie  original  values  of   r 
and  r,      Ax-j  (cp)  being  used  instead  of   ZiX^(cp),  the  re- 
sulting chordv/ise  extremities  would  in  general  not  be  at 


X  =   ±1.   It  is  therefore  necessary  to  adjust   r 
such  that  v;ith  the  derived  L\x^,    Ay]_, 


and 


X- 


X- 


-r%)  =  ' 


=  -1 


(21^) 


where  <^j        and   c,o„   are  the  axigles  on  the  circle  corre- 
sponding to  the  extremities  of  the  desired  airfoil.   This 
operation  was  mentioned  in  the  section  "Superposition  of 
Solutions."   It  may  be  termed  a  horizontal  stretching  of 
the  given  airfoil.   The  condition  given  by  equations  .  (.2ii.) 
applied  to  equation  (9)  yields 


1  =  T^  +  r3_  cos  (Pj.  +  Ax]_^cp 


Ni) 


> 


-1 


T3_  +  ri   cos  <?rp   +  AXt /Q^  "N 


(25) 


Subtraction  of  the  second  of  these  equations  from  the 
first  gives  for   r-. 


1  + 


^^^i(^t-l)  -  ^^i(y^) 


cos  cp- 


Ni 


-  cos  cp, 


(26) 


Tl 
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Addition  of  equations  (25)  gives  for   T, 

cos  cpTj   +  cos  cp^    '^"f'^N  )  "^  ^^'^if'^T   ) 

• +   -S : .; i_ 


"l   = 


(2?) 


The  angles  (D-^        and  Om    in  equations  (26)  and  (27) 

correspond  to  the  extrerriities  of  the  desired  airfoil. 
They  are  given  by  graphical  solution  of  equation  (11) 


dZix-i  (cp) 

sin   cp  = 

r2_   dcp 


(11) 


Equation  (11)  must  be  solved  simultaneously  with  equa- 
tion (26)  for   r]_,   ""-Pit  J   and 

few  successive  trials  ere  necessary.   Thence   t 


cpm,  .   In  practice  only  a 


IS 


obtained  by  equation  (27)  .   The  angle   cPrn    determined  in 

this  process  is  equivalent  to  the  zero-lift  angle  of  the 
airfoil,  equation  (12). 


Illustrative  Example  of  Direct  Method 


As  a  numerical  Illustration  of 
velocity  distribution  of  the  NAGA  65 
culated.  In  order  to  obtain  an  inlt 
of  the  6-percent  camber  circular  arc 
was  added  to  the  CMP  of  the  12-Derce 
profile,  derived  from  that  of  table 
previous  section.  Before  this  addit 
CLIP  for  the  circular  arc  v/as  Increas 
plied)  by  1.0928/1 .OO72  to  correspo 
of  straight  line  r  as  the  sjaimetri 
normalized  resultant  CI\IP  and  the  ass 
given  in  tables  V(a)  and  VI,  respect 
airfoil  is  shov/n  in  figure  7(a.). 


the  direct  method  the 
12  airfoil  v/as  cal- 
lal  airfoil,  the  CMP 

(tables  III  and  IV) 
nt  thick  syi-nmetrical 
I  as  indicated  in  a 
ion  v/as  made ,  the 
ed  in  scale  (multi- 
nd  to  the  same  length 
cal  profile  ClIF .      The 
ociated  constants  are 

vely.   The  initial 


The  given  airfoil,  IJ3CA  6512,  was  so  rotated  through 
an  angle  of  -0.88"^  (nose  dov/n)  as  to  be  tangent  to  the 
initial  airfoil  at  the  leading  edge.   The  convergence 
near  the  leading  edge  v/as  thereby  accelerated.   The  given 
airfoil  is  shov/n  in  this  position  in  figure  7(a-)«   Two 
approximations  v/ere  then  carried  out  in  accordance  with 
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the  procedure  given  in  the  preceding  section.   The  nuineri- 
cal  results  are  given  in  tables  V  and  VI.   The  first 
approximate  airfoil  is  indicated  by  the  circles  in  fig- 
ure 7(^)j  the  second  approximate  airfoil  was  indistin- 
guishable to  the  scale  used  (chord  =  20  in.)  fron  the 
given  airfoil.   The  velocity  distributions  of  the  initial, 
first,  and  second  approximate  airfoils  are  given  in  fig- 
ure 7(^)>  together  with  those  corresponding  to  one 
approximation  by  the  Thecdorsen-Garrick  method  (refer- 
ence 5)«   J^hs  second  approximation  velocity  distribution 
differs  appreciably  from  that  of  the  Theodorsen-G-arrick 
method  on  the  UDper  surface  but  agrees  faiz-^lv  well  on  the 
lov/er  surface.   The  discrepancy  for  the  rearmost  S  percent 
of  chord  on  the  lower  surface  appears  to  be  due  to  lack 
of  detail  in  this  region  in  the  Theodorsen-Garrick  cal- 
culation. 


The  convergence  of  the  CMP  method  is  seen  to  be 
rapid,  considering  the  approximate  nature  of  the  initial 
airfoil,  although  two  approxim.ations  are  required  for  a 
satisfactory  result.   The  second  approximation  could 
probably  have  been  made  unnecessary  by  suitably  adjusting 
the  first  increment   5y-|  (cp)   near  the  leading  and  trailing 

edges  on  the  upper  surface   before  calculating   5xt_(CD). 

The  direction  in  v/hich  to  adjust  the  increm.ent  is  obtained 
by  comparing  the  thicimess  of  the  initial  airfoil  v:ith 
that  of  the  given  airfoil  in  these  regions.   Because  a 
thicker  section  has  a  greater  concentration  of  ohordv/ise 
locations  tcv/ard  the  extremities,  for  a  given  set  of 
cp   points,  than  does  a  thinner  section,  the  chordwise 
stations  woiald  be  expected  to  be  shifted  outward  as  the 
thickness  of  the  section  is  increased.   The  ordinates 
'i^yi(^-P)   should  therefore  have  been  chosen  at  chordwise 
stations  slightly  more  toward  the  extreir^ities  than  those 
given  hj   equation  (9)- 

The  accuracy  of  the  velocities  is  estimiated  to  be 
within  1  percent.   It  was  expected,  and  verified  by  pre- 
liminary calculations,  that  the  results  v/ould  tend  to 
be  more  Inaccurate  tov/ard  the  extremities  of  the  airfoil 
than  near  the  center.   This  result  is  evident  fromi  equa- 
tion (17).   A  given  inaccuracy  in  the  slopes   dax/d^c   and 
dAy/dcp   can  produce  a  large  error  in  the  velocity  near 
the  extremities,  v/here   sin  Cp   approaches  zero.    This 
disadvantage  does  not  appear  in  the  Theodorsen-Garrick 
method,  in  v.fhich   sin  cp-   is  replaced  by  one.    Excessive 
error  in  these  regions  can  be  avoided  in  various  ways. 
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If  the  initial  airfoil,  for  which  the  slopes   dAx^/dCp 
and   dAyo/di?  have  presumably  been  computed  accurately, 
Is  a  good  approximation  in  these  regions,  as  evidenced 
by  the  smallness  of   6X]_,   5y]_   compared  to   Ax^,   Ayo> 
the  effect  of  Inaccuracy  of  the  slopes   dox^L/'^^^j  d5y-j_/dcp 

will  be  reduced,  since  they  are  added  to  the  initial 
slopes   dAXp/dT',   dAy^/dcp.   It  was  to  reduce  the  magnitude 
of  the  incremental  CI'.IP  near  the  leading  edge  that  the 
IIACA  6512  airfoil  was  drawn  tangent  to  the  initial  air- 
foil in  this  region. 

The  error  in  the  derivatives  can  also  be  avoided  by 
computing  them  from  the  differentiated  Fourier  series 
for   5x-]^,   Gy]  .   (See  appendix  A.)   This  calculation  was 
made  in  the  illustrative  example,  after  it  was  found  that 
an  error  of  about  5  psrcent  in  tlie  velocity  on  the  upper 
surface  leading  edge  could  be  caused  by  unavoidable 
inacciiracy  in  measuring  the  Incremental  slopes. 

The  fact  that  the  comptited  derivatives  do  not  repre- 
sent the  derivatives  of  the  CMF  but  rather  the  deriva- 
tive of  its  Fourier  expansion  to  a  finite  niurriber  of 
terms  may  introduce  inaccuracy.   (The  derivative  Fourier 
series  converges  more  slovvly  than  the  original  series.) 
A  com.parison  of  the  computed  derivatives  vn.th  the  measured 
slopes  will  Indicate  the  limits  of  error,  however,  as  v.'ell 
as  the  true  derivative  curve. 

The  importance  of  knov/lng  the  CMP  derivatives  ac- 
curately may  make  it  desirable  to  solve  the  direct 
problem  from  the  airfoil  slopes,  rather  than  from  the 
airfoil  itself,  as  given  data.    This  variation  of 
technique  enables  the  CMF  derivatives  rather  than  the 
ClIF   itself  to  be  approximated  initially.    Further 
details  are  given  in  reference  1. 


II  -  THE  INVERSE  POTENTIAL  PROBLEM  OF  AIRFOIL  THEORY 


The  inverse  potential  problem,  of  airfoil  theory  may 
be  stated  as  follows:   Given  the  velocity  distribution 
as  a  function  of  percent  chord  or  surface  arc  of  an  unknown 
airfoil  -  to  derive  the  airfoil.   Before  the  questions  of 
existence  and  uniqueness  of  a  solution  to  the  problem  as 
thus  stated  are  discussed,  several  CMP  methods  of  solu- 
tion will  be  outlined  and  illustrated  by  numerical 
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examples.   Various  previous  methods  of  solution  will  then 
be  described  briefly  and  their  inherent  limitations  and 
restrictions  on  the  prescribed  velocity  distribution  will 
be  compared  with  those  of  the  CMF  m.ethods. 

The  prescribed  velocity  distribution  is  assiomed  to 
be  either  a  double-valued  continuous  function  of  the 
percent  chord  or  a  single -valued  continuou.s  function  of 
percent  arc.   (Isolated  discontinuities  in  velocity  are, 
hov:/ever,  at  least  in  the  percent-chord  case,  admissible.) 


CMF'  Method  of  Potentials 

This  Inverse  method  is  based  on  the  fact  that,  if 
the  airfoil  and  its  corresponding  flat  plate  and  circle 
are  immersed  in  the  sane  free -stream  flows  and  have  the 
same  circulation,  conformally  corresponding  points  in 
the  three  planes  have  the  saine  potential. 

Consider  first  the  case  where  a  velocity  distribu- 
tion corresponding  to  a  s7,Tm-fietrloal  airfoil  at  zero  lift 
is  specified  as  a  function  of  percent  chord.   If  an 
initial  airfoil  is  assumed,  the  prescribed  velocity  can 
be  integrated  along  its  surface  to  3^ield  an  approximate 
potential  distribution  as  a  function  of  percent  chord. 
This  potential  increases  from,  zero  at  the  leading  edge 
to  a  maxira,um  value  at  the  trailing  edge.   Of  fundamental 
importance  to  the  success  of  the  method  is  the  fact  that 
this  potential  cm'^ve  depends  mainly  on  the  prescribed 
velocity  distribution  and  only  to  a  much  lesser  extent 
on  the  form  of  the  initially  assumed  airfoil.   The  chord 
line  of  the  initial  sirfoil  taken  as  the  x-axis  is  next 
sufficiently  extended  that,  in  the  same  free-stream  flow 
as  for  the  airfoil,  the  potential,  which  in  this  case 
is  sim.ply   V£,  increases  linearly  from  zero  at  its 
leading  edge  to  the  same  miaximtim  value  at  the  trailing 
edge  as  exists  for  the  approximate  potential  curve  derived 
initially.   Horizontal  displacements   ax   betv-/een  these 
curves  are  then  measured  as  a  function  of  the  straight- 
line  abscissas  and,  hence,  as  a  function  of  the  central 
angle   cp  of  ■  the  circle  corresponding  to  the  straight 
line.   These  horizontal  displacements   /ix(cp),  together 
with  the  conjugate  function  Ay(cp)   computed  therefrom 
and  the  length  of  straight  line  previously  determined, 
constitute  a  CMP  for  an  airfoil  that  is  a' first  approxi- 
m.ation  to  the  unknown  airfoil.   The  approximation  is 
based  on  the  use  of  a  more  or  less  arbitrarv  initial 
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airfoil  to  set  up  the  first  approximate  potential.   The 
exact  velocity  distribution  of  the  derived  fii'-st  approxi- 
mate airfoil  can  now  be  coraputed  and  compared  vn*.th  the 
prescribed  velocity.   If  the  agreement  is  not  satisfac- 
torily close,  the  procedure  is  repeated,  with  the  airfoil 
just  derived  taking  the  place  of  the  one  initially  assumed. 

The  complication  introduced  in  the  general  case  in 
which  the  prescribed  velocity  distribution  corresponds 
to  an  ursyimrietrlcal  airfoil  with  circulation  can  be 
resolved  as  follows:   It  is  convenient  in  this  case  to 
discuss  the  potentials  in  the  circle  plane.   The  pre- 
scribed velocity  distribution  is  transferred  to  the  circle 
plane  by  means  of  the  stretching  factor,  presiomed  known, 
of  the  initially  assumed  airfoil;  that  is,  equation  (li;) 
is  solved  for  v,-^((P).   The  first  approximate  potential 
distribution  as  a  function  of  the  central  angle   cp   is 
obtained  by  integrating  Vq(o)   through  £i  <^P-range  of  2.rr 

radians  (around  the  airfoil),  starting  from  the  value 
of  <?■   near  zero  for  which  Vp(<P)   is  zero  (the  front 

stagnation  point).   This  approximate  potential  curve  has 
a  minlmiua  value  of  zero  at  the  front  stagnation  point, 
rises  to  a  maximum  for  the  value  of   cp  near   rr   corre- 
sponding to  the  rear  stagnation  point,  then  falls  to  a 
minimuTfi  for  the  final  value  of   cp   (the  front  stagnation 
point),  which  is  an  angle   2tt  radians  from  the  starting 
cp-point.   The  difference  between  the  final  and  the  initial 
potential  minimums  is  a  first  approximation  to  the  circu- 
lation r. 

A  circle  of  such  diameter  is  now  derived  v/hich,  with 
this  circulation  and  the  same  free-stream  flow  as  for  the 
airfoil,  yields  a  |)otential  distribution  (henceforth  called 
true  potential  distribution)  that  has  the  same  maxim.um 
and  minimiJim  values  as  the  approximate  potential  curve 
just  derived.   If  the  maximum  approximate  potential  is 
denoted  by  rQU  and  the  decrease  of  potential  (considered 
positive)  from  the  maximurii  to  the  final  value  by  v^L, 
where  Vq      is  the  diameter  of  the  circle  corresponding 
to  the  initial  airfoil,  the  parar:ieter   y  i^  first  com- 
puted from 


TT        U  -  L 


2(y  +  cot  y)    U  +  L 


(28) 
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by  means  of  figure  8.   The  desired  diameter   r   is  then 
given  by 

ro(U  +  L) 

r  =  7^ : 7  (29) 

1^(  cos  Y  "^  Y  sm  y) 

The  parameter   y   i^  actually  the  sum  of  the  angle  of 
attack  and  zero-lift  angle  of  the  unknown  airfoil,  to  a 
first  approximation;  that  is, 

Y  =  a  +  3^  (50) 

It  is  related  to  tlie  circulation  T      by  equation  (l8). 

This  procedure  for  the  calculation  of  the  diameter 
(see,  for  example,  reference  6)  follows  easily  from,  the 
expression  for  the  poteiitial  distribution  on   a  circle, 
obtained  by  integration  of  equation  (15)  ^s, 

$t(®)  =  i-n  /   v^((;o)dCp 


t- 


^o 


[cos  Y  "'"  Y  sin  y  ~  ^c>s  (<y'^-  a)  +  ( c?+  a)  sin  y|    (31) 


If  the  diajiieter   r   of  the  derived  circle  is  much 
greater  than  the  diameter   r^   of  the  013.-^010  corresponding 
to  the  initial  airfoil,  it  is  desirable  to  increase  the 
CMP  AXq,   Ay^   of  the  initial  airfoil  by  a  factor  suffi- 
cient to  m.odify  the  initial  airfoil  such  that  it  corre- 
sponds to  a  circle  of  diameter   r.   A  new  approxim.ate  and 
true  potential  distribution  is  then  obtained  as  described 
but  by  using  the  modified  initial  airfoil. 

The  first  apprcxirmte  horizontal  displacem.ent  func- 
tion is  now  determined  as  the  sum  of  the  horizontal 
displacement   AXq(cp)   corresponding  to  the  (raodified) 
Initial  girfoll  and  an  Increment   Sx-,  (O))   produced  by 
the  noncoincidence  of  th.e  approximate  potential  distri- 
bution ^g^   and  the  true  potential  distribution  ^^. 

This  horizontal  increment  may  be  m.easured  betv/een  the 
two  potential  curves,  both  considered  plotted  against 
chordwlse  position  in  the  physical  plane.   With  sufficient 
accuracy  this  increment  may  be  computed  as  the  vertical 
distance  betv/een  the  potential  curves  divided  by  the 
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airfoil  to  set  up  the  first  approximate  potential.   The 
exact  velocity  distribution  of  the  derived  first  approxi- 
mate airfoil  can  now  be  computed  and  compared  with  the 
prescribed  velocity.   If  the  agreement  is  not  satisfac- 
torily close,  the  procedure  is  repeated,  with  the  airfoil 
just  derived  taking  the  place  of  the  one  initially  assumed. 

The  complication  introduced  in  the  general  case  in 
which  the  prescribed  velocity  distribution  corresponds 
to  an  unspnmetrical  airfoil  with  circulation  can  be 
resolved  as  follows:   It  is  convenient  in  this  case  to 
discuss  the  potentials  in  the  circle  plane.   The  pre- 
scribed velocity  distribution  is  transferred  to  the  circle 
plane  by  means  of  the  stretching  factor,  presumed  knov/n, 
of  the  initially  assuraed  airfoil;  that  is,  equation  {l\.\.) 
is  solved  for  v^(c,o).   The  first  approximate  potential 
distribution  as  a  function  of  the  central  angle   Q   is 
obtained  by  integrating  VQ(cr')   through  a  cp_range  of  Zrr 

radians  (around  the  airfoil),  starting  fx'^om  the  value 
of  <¥   near  zero  for  which  Vp(cp)   is  zero  (the  front 

stagnation  point).   This  approximate  potential  curve  has 
a  mlnlmujn  value  of  zero  at  the  front  stagnation  point, 
rises  to  a  maximiun  for  the  value  of   c?  near  v     corre- 
sponding to  the  rear  stagnation  point,  then  falls  to  a 
minimum  for  the  final  v?lue  of   cp   (the  front  stagnation 
point),  v.'hich  is  an  angle   2rr  radians  from  the  starting 
cp-point.   The  difference  between  the  final  and  the  initial 
potential  minimums  is  a  first  approximation  to  the  circu- 
lation r. 

A  circle  of   such  diameter  is  now  derived  which,  with 
this  circulation  and  the  same  free-stream  flov/  as  for  the 
airfoil,  yields  a  potential  distribution  (henceforth  called 
true  potential  distribution)  that  has  the  same  maxim.um 
and  minimijun  values  as  the  approximate  potential  curve 
just  derived.   If  the  maximum  approximate  potential  is 
denoted  by  rQU   and  the  decrease  of  potential  (considered 
positive)  from  the  maximum  to  the  final  value  by  VqL, 
v/here   ro   is  the  diameter  of  the  circle  corresponding 
to  the  initial  airfoil,  the  parameter   y  is  first  com- 
puted from 


TT         U  -  L 


2(y  +  cot  y)    U  +  L 


(28) 
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by  means    of   figuz'-e    8.      The    desired   diaraeter      r      is    then 

.    To^U  +   L) 


given    uy 


Li_(  cos    Y  "^    Y    sin    y) 


(29) 


The   parameter      y      ^^    actually   the    sura   of    the    angle    of 
attack  and   zero-lift   angle   of   the   -unknov-fn   airfoil^    to   a 
first   approximation;    tiiat   is, 

r  =   a  +   Q^  (50) 

It  is  related  to  the  circulation  F   by  equation  (l8). 

This  procedure  for  the  calculation  of  the  diameter 

(see,  for  example,  reference  6)  follows  easily  from  the 

expression  for  the  potential  distribution  on  a  circle, 
obtained  by  integration  of  equation  (15)  as 


o  /   v,^(co)d?) 


^o 


fcos  Y  +  Y  sin  y  -  '^'^^   ('*' ^  a)  +  (c?+  a)sin  y|    (51) 


If  the  diajneter   r   of  the  derived  circle  is  much 
greater  than  the  diameter  r^   of  the  circle  corresponding 
to  the  initial  airfoil,  it  is  desirable  to  increase  the 
CMP  AXq,   Ay(3   of  the  Initial  airfoil  by  a  factor  suffi- 
cient to  modify  the  initial  airfoil  such  that  it  corre- 
sponds to  a  circle  of  diameter   r.   A  nev/  approximate  and 
true  potential  distribution  is  then  obtained  as  described 
but  by  using  the  modified  initial  airfoil. 

The  first  approximate  horizontal  displacement  func- 
tion is  nov/  determined  as  the  sum  of  the  horizontal 
displacement  AXg(cp)   corresponding  to  the  (modified) 
initial  airfoil  and  an  increment   Ox-,  (®)   produced  by 
the  noncoincidence  of  the  approximate  potential  distri- 
bution ^g^      and  the  true  potential  distribution  0^. 

This  horizontal  increment  may  be  measured  betvifeen  the 
two  potential  curves,  both  considered  plotted  against 
chordvvise  position  in  the  physical  plane.   V/ith  sufficient 
accuracy  this  increment  m^ay  be  computed  as  the  vertical 
distance  between  the  potential  curves  divided  by  the 
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slope  of  the  spproxiinate  potential  curve;  namely,  the 
prescribed  velocity  v  .   If,  thei 
are  considered  as  functions  of   cp 


prescribed  velocity  v  .   If,  therefore,  all  quantities 


Ax^^  =  AXq  +  5x-[_ 

0  (CD)  -  0  (cp) 
=   AXq  +  -^ — ^ 52 

The  ordinate  function  Ay-i(C)   conjugate  to   Ax-il"?) 

can  nov/  be  coniputed  and,  together  with   Ax-l(cp)   and  the 
diameter   r   obtained  previously,  determines  the  first 
approximate  airfoil  by  equations  (9)  and  (10).    Calcu- 
lation or  measurement  of  the  CMP  derivatives   dAx2_/dcp, 
dAy^^/d©  and  the  use  of  equations  (11)  and  (17)  then 
determine   the  zero  lift  angle   '3    and  the  exact  velocity 

distribution  of  the  first  ajjproximate  airfoil.   The  angle 
of  attack,  to  a  first  approximation,  is  given  by  equa- 
tion (50) J  the  value  of   y  derived  from  equation  (23) 
being  used.   This  exact  velocity  distribution  is  compared 
with  that  prescribed  and,  if  the  agreement  is  not  close 
enough,  the  procedure  can  be  repeated  with  the  first 
approximate  eiirfoil  as  the  initial  airfoil. 

In  the  case  where  the  prescribed  velocity  is  speci- 
fied as  a  function  of  percent  arc,  then  by  line  integra- 
tion of  the  prescribed  velocity  along  the  percent  arc, 
the  true  potential  distribution  of'  the  unknown  airfoil 
is  known  as  a  function  of  arc  (except  for  a  trivial  scale 
factor).   The  maximum  and  ninimum  values  of  this  potential 
distribiition  then  perm.it  the  'unique  determination,  by  the 
calculation  T^rsviously  described,  of  the  circle  corre- 
sponding conformally  to  the  unknovvn  airfoil.   Correlation 
of  the  potential  distribution  of  this  circle  v/ith  the 
potential  distribution  as  a  function  of  arc  initially 
calculated  therefore  7/ields  exactly  the  potential  distri- 
bution of  the  unknovm  airfoil  as  a  function  of  the  central 
angle  c?   of  the  circle.   This  fact  has  been  noted  by 
Gebelein  (reference  6),  The  calculation  of  the  diameter   r 
as  outlined  above  for  the  percent-chord  case  is  thus 
unnecessary.   The  remainder  of  the  procedure  is  the  same, 
the  successive  approximate  airfoils  now  being  adjusted 
to  correspond  conformally  to  this  circle,  before  corre- 
lating their  percent -arc  lengths  v»lth  the  prescribed 
velocity  distribution  in  preparation  for  the  next 
approximation. 
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The  successive  contours  determined  by  the  method  of 
potentials  are,  of  necessity,  closed  contours,  whether  or 
not  the  sequence  of  contours  converges  to  a  solution 
satisfying  (mather.iatically)  the  prescribed  velocity  dis- 
tribution.  The  closure  of  the  contours  is  a  consequence 
of  the  method  of  setting  up  the  horizontal  displace- 
ments,  A>:(p),   and  solving  for   Ay(cp),  by  wMch  the 
contour  coordinates  are  obtained  as  single-valued  func- 
tions of   <9.   The  necessity  for  closed  contours  does  not, 
however,  exclude  the  possibility  of  deriving  physically 
unreal  shapes;  namely,  contours  of  figure-eight  type. 
This  point  will  be  discussed  at  greater  length  later  but 
it  may  be  remarked  here  that  it  is  the  extra  degree  of 
freedom  introduced  by  the  class  of  figure-eight  type 
contours  that  adinits  the  possibility  of  a  unique  solu- 
tion to  the  inverse  problem  treated  here. 

It  will  have  been  noticed  that,  whereas  in  the  direct 
method  a  l\j      is  determined  from  the  given  data  -  that  is, 
the  airfoil  -  and  a   Ax   is  computed  therefroFi,  conversely, 
in  the  inverse  method  of  potentials  a   Ax   is  determined 
from,  the  given  data  -  that  is,  the  velocity  distribution  - 
and  a   Ay   is  computed  therefrom.   Slm.ilarly,  just  as 
the  direct  problem  can  also  be  solved  by  deriving   dAy/dCp 
from  the  given  airfoil  slopes  and  thence  computing 
dAx/dcp,  so,  conA-ersely,  can  the  Inverse  problem  be  solved 
by  deriving   dAx/dcp   from  the  prescribed  velocity  dis- 
tribution and  thence  computing   dAy/dcp.   This  inverse 
method  of  derivatives  will  be  discussed  after  som.e 
numerical  examples  are  presented.  Illustrating  the  m.ethod 
of  potentials. 


Exam.ples  of  GI.TP  Method  of  Potentials 

Symmetrical  section.  -  The  m.ethod  of  potentials  was 
applied  first  to~tKe  derivation  of  the  s^/mmetrical  profile 
corresponding  to  the  prescribed  velocity  distribution 
shown  in  figure  9(a).   As  an  initial  airfoil  the  12- 
percent  thick  profile  derived  from  the  5'3-percent  thick 
Joukowskl  profile  in  part  I  was  used.   The  initial  CLIF 
and  associated  constants  are  given  in  table  VII.   The 
initial  airfoil  and  its  velocity  distribution  are  shov/n 
in  figure  9.   The  first  increment  CMP  and  the  resultant 
first  approximate  airfoil  and  its  exact  velocity  distri- 
bution v^rere  calculated  by  the  procedure  of  the  preceding 
section.   The  increm.ental  slopes   d5x]_/d(p,   d5y]_/dcp 
were  computed  and  found  to  approximate  the  measured  slopes 
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very  closely.   The  results  are  presented  in  table  VII 
and  figure  9.      It  is  seen  that  the  change  in  velocity 
and  profile  accomplished  by  one  step  of  the  Inverse 
process  is  large;  that  is,  the  convergence  is  rapid. 
The  high  velocity  of  the  first  point  on  the  upper  surface 
(cp  =  15°)  is  due  to  lack  of  detail  in  the  calculation. 
(Twelve  points  on   the  upper  surface  were  calculated.) 
Por  practical  purposes  the  nose  could  be  easily  modified 
to  reduce  this  velocity  if  desired  without  going  through 
a  complete  second  appro::imatlon. 

Mean  camber  line  for  uniform  velocity  increiaent . - 


As  a  second  example  of  the  inverse  CMF  method,  the  profile 
producing  uniform  equal  and  opposite  velocity  increments 
on  upper  and  lov;er  surfaces  was  derived.   By  the  methods 
of  thin-airfoil  theory  this  velocity  distribution  yields 
the  so-called  logarithmic  camber  line.   The  prescribed 
velocity  distribution  is  indicated  in  figure  10(a).   The 
velocity  peaks  at  the  extremities  of  the  prescribed 
velocity  curve  ivere  assumed  in  order  to  compensate  for 
an  expected  rounding  off  of  the  velocity  in  this  region 
in  working  up  from  the  initial  velocity  distribution. 
The  convergence  to  the  prescribed  uniform  velocity  dis- 
tribution would  thereby  be  accelerated.   The  initial 
airfoil  was  taken  as  the  o-percent  camber  circular  arc, 
discussed  in  part  I.   The  initial  CIvIP  and  its  associated 
constants  are  given  in  tables  III  and  IV.   The  circular 
arc  and  its  velocity  distribution  are  shov/n  in  figure  10. 


A  first  approximation  v/as  calculated  as  outlined  in 
the  previous  section.   A  niAi.ierical  difficulty  appeared 
in  the  process  of  solving  equation  (11)  for  the  zero- 
lift  angle  of  the  first  approximate  airfoil.   It  appeared 
that  a  2k-polnt  calculation  (12  points  by  sjTmetry)  did 
not  give  sufficient  detail  in  the  range   tt  <  cp  <  13.  tr 

to  yield  a  reliable  solution  of  equation  (11)  for  the 
zero-lift  angle.   This  result  was  a  consequence  of  the 
prescribed  velocity  discontinuity  et  the  extremities  with 
the  consequent  large  but  local  changes  in  CliP  and  profile 
shape  required  in  these  regions.   The  solution  obtained 
for  the  zero-lift  angle  was   Pm  =  6.1°,   v/hich  by  equa- 
tion (19)  with   r  =  I.OOJ4.5   and   a3_  =  0   yielded 
cj  =  0.67.   The  desired   c,,   hov.'ever,  is  O.8O,  which 
would  correspond  to  ^^  =   7.27°.   It  was  considered 

that  a  relatively  minute  change  in  the  shape  of  the 
extremities  of  the  derived  camber  line  vvould  alter  the 
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slope   dAx-[_/dcp   In  the  desired  range  sufficiently  to 
yield  a  zero-lift  angle  of   Pm  =  7.27°.   On  the  other 

hand  the  effect  of  such  a  local  change  on  the  CMF  as  a 
whole  would  be  snail.   The  velocity  distributions  of  the 
derived  profile  v/ere  therefore  computed  for  both  zero- 
lift  angles  quoted  previously. 

The  results  are  given  in  table  VIII  and  in  figure  10. 
Included  for  comparison  in  figure  10(b)  (vertical  scale 
magnified)  is  the  logarithmic  mean  line  of  thin-airfoll 
theory,  computed  for  c-,    -   0.80.   The  velocity  distri- 
bution of  the  derived  shape  as  calculated  for  the  desired 
lift  coefficient  of  c-,    -   0.80   is  seen  to  be  a  satis- 
factory approximation  to  the  desired  rectangular  velocity 
distribution.   The  profile  itself  is  seen  to  be  one  of 
finite  thickness  as  compared  with  the  single  line  of  ■ 
thin-airfoil  theory. ■  Airfoils  obtained  by  superposition 
of  this  type  of  camber  line  with  thickness  profiles  would 
therefore  he    increased  in  thickness  over  that  of  the 
basic  thickness  form. 

The  changes  in  velocity'-  distribution  and  in  shape 
of  profile  are  again  seen  to  be  large j  that  is,  the  con- 
vergence was  rapid.   As  is  to  be  expected,  the  rapidity 
of  convergence  of  both  the  direct  and  inverse  methods  in 
comparable  cases  is  about  the  same. 


CMF  Method  of  Derivatives 

Instead  of  approximating  by  the  method  of  potentials 
to  a  CMF  that,  when  differentiated,  yields  the  prescribed 
velocity,  the  CMP  derivatives  may  be  obtained  directly. 
The  controlling  equations  are  equations  (17)j  (9)j  s^^d 
a  modification  of  equation  (7)« 

I  s In    ( CO   +    a )    +    s  in    (a  +    Pm)  i 
vz(cp)    =    ' _  ^         n^  (17) 

\ri^    .    3ln   cpf   .    (-^f 
V\r   dcp  /  Vr  dcp/ 


nSrr 

dAy  _     1      I  dAx         ,    <P'    -    cp 


dcp  2Tr  dcp' 

Jo 


cot   ^ ^    dcp'  (7a) 
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^  =   cos    cp     +  -  Ax^o)  (9) 

r  r 

These  equations,  together  v/ith  the  aui^iliary  equations  (11) 
and  (l8),  constitute  a  set  of  s  inailtaneous  equations  from 
Vvhlch  the  CIIF  dei-^lvative   dAx/dc?  may  be  determined  from 
a  prescribed  volcclt;'-  distribution   v^ .   The  corresponding 
airfoil  is  determined  by  integration  of   dAx/dcp  and  its 
conjugate   dAy/dcp. 

Consider  first  tbe  case  where  bhe  velocity  is  speci- 
fied as  a  function  of  percent  arc.   As  explained  in  the 
previous  section,  the  constants   r   and   y   of  "the  final 
circle  corresponding  to  the  unknovm  airfoil  can  in  this 
case  be  determined  initially.   Points  of  equal  potential 
along  the  arc  and  circle  are  then  found,  which  yield   Vg 
as  a  function  of   cp.   The  angle  of  attack   a   in  equa- 
tion (17)  Is  taken  as  sor.ie  reasonable  value  and   dAx/r  dcp 
determined  by  successive  approximation.   In  the  first 
approximation   dAy/r  dcp  nay,  for  example,  correspond  to 
some  known  Gll^ .      Equation  ( I7 )  is  then  solved  for 
dAx/'r  dcp,   for  v/hich  the  conjugate   dAy/r  clcp  is  calcu- 
lated next  and  used  as  a  basis  for  a  better  determination 
of   dAx/r  dcp.   The  airfoil  corresponding  to  any  approxi- 
mation Is  obtained  by  integration  of   dAx/dCf'  and  its 
conjugate   dA^v'^^cp.   (The  method  of  derivatives  may  be 
regarded  as  based  on  the  use  of  the  function 

i-D  — — .   This  function  is  regular  everywhere  outside 

ICP 

the  circle   p  -  Re   ,   approaches  zero  at  infinity,  and 

,    dAx  .  .  dAv      ,  ,       -I   .  ,  T  ^     \ 
reauces  to  -- —  +  1  —7-^  on  tJie  circle  itself.  ) 
Q  cp      dco 

In  general  the   dAx/dcp  as  determined  in  any  approxi- 
Fiatlon  will  have  an  average  value  other  than  aero.   The 
Ax(ci:')   obtained,  say,  by  integration  of  its  Fourier 
series  would  therefore  contain  a  term  proportional  to  co 
in  addition  to  a  Fourier  sez-'ies.   Thus,   Ax(cp)   would 
not  be  a  single-valued  function  of   cp   and  the  resulting 
contour  would  not  close.   Simply  subtracting  the  average 
value  of   dAx/c3cp   (the  constant  term  in  Its  Foiirier  series), 
however J  will  close  the  derived  contour.   If  the  method 
converges,  this  average  value  approaches  zero  in  the  suc- 
cessive approximiations . 

A  preliminary  over-all  adjustment  of  an  initially 
chosen  CIvIF  ma;/  be  desirable.   Thus,  if   dAx-j_/dcp   is 
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calculated  in  terr,is  of  the   dAyQ/do?   of  a  previous  approxi- 
mation and  is  found  to  be  larger  than   dAx^/dcp   by  some 
factor,   dAyg/dcp   can  be  multiplied  by  this  factor  and 
the  calculation  of  dl\x-^/d(p      repeated. 

Although  the  angle  of  attack  may  be  arbitrarily  set 
initially  in  this  calculation  it  should  be  so  chosen  that 
the  final  airfoil  Virill  coincide  approximately  in  position 
v/ith  the  initial  airfoil.   After  each  calculation  of 
dAx/dcp,   the  zero-lift  angle   Pm   can  be  calculated, 

equation  (11),  which  thereupon  fixes   a,   since   y  -  °- "*"  Pm 
is  known. 

If  the  prescribed  velocity  distribution  is  speci- 
fied as  a  function  of  percent  chord,   v^lcr)   must  be 
determined  in  the  successive  approximations  by  use  of 
equation  (9)'   The  quantity  f^-   =    a  +    ^^^     may  be  deter- 
mined in  each  approximation  as  in  the  method  of  potentials 
or,  in  physically  real  cases,  by  equation  (19).   The 
diameter   r   is  so  determined  that  the  successive  airfoils 
are  of  a  standard  chord  length. 

It  is  evident  from  the  structure  of  equation  (17) 
that  near  the  airfoil  extremities  v;here   sin  cp — >  0,  and 
in  particular  at  the  nose  of  the  airfoil  v/here   dAy/dcp 
is  comparable  to   dAx/dcp   in  m^agnitude,  the  convergence 
by  this  method  (and  by  the  method  of  potentials)  will  be 
comparatively  slow.   If  modifications  to  the  airfoil  only 
in  the  immediate  neighborhood  of  the  nose  are  required, 
it  may  be  more  expedient  to  apply  a  preliminary  Joukowski 
transformation,  that  is,  to  use  these  methods  with  the 
Theodorsen-Garrick  transf ormiation. 

An  example  of  the  use  of  the  ClIF   method  of  deriva- 
tives to  solve  an  inverse  problem  is  given  in  reference  1 
for  the  case  of  a  cascade  of  airfoils. 


Method  of  Betz 

In  the  inverse  method  of  Betz  .  (reference  7)  slu   air- 
foil and  its  velocity  distribution  are  assumed  known 
(fig.  11)  and  a  desired-  velocity  is  specified  as  a  func- 
tion of  percent  arc.   The  nev/  velocity  and  length  of  arc 
are  specified  in  such  a  way  that  the  extremities  of 
potential  are  the  sam.e  as  on  the  known  airfoil.  .  Both 
known  and  unknown  airfoils  then  transform  into  the  same 
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circle  and,  in  particular,  the  velocities  at  points  of 
equal  potential  on  the  two  profiles  can  be  found. 

In  order  to  deterruine  the  profile  corresponding  to 
the  new  VGlccity,  the  complex  displacement   Zo  -  z-, 
between  points  of  equal  potential  on  the  two  profiles  is 
expressed  as  a  function  of  the  CGrr3sponding  complex 
velocities  (denoted  07  Vg )  thus. 


dz- 


I J  dZj_ 


d^.'/dz^        ^Zj_ 


aw 


'7dz 


2 


Vr. 


Hence 


Zo  -  z 


.33) 


where  the  integration  is  cerrled  out  elong  the  known  pro- 
file from  the  trailing  edge,  which  is  taken  es  coincident 
for  the  two  airfoils,  to  the  point   z-i.   The  complex 


I'unction 


V, 


1/^-2 


known  rstio 


V. 


is  determined  approximately  from  the 


corresponding  to  the  points  of  equal 


potential  by  the  argument  that,  inasmuch  as  the  two  pro- 
files hsve   nearly  the  s sme  slope  at  corresponding  points. 


the  real  iDart  of 


V. 


V, 


n 


1   is  given  by 


■1 


v. 


-  1.   (This 


assumption,  like  the  approximations  in  the  CMP  methods, 
is  least  valid  at  the  nose  of  the  airfoil.   The  function 
Zp  -  z-|   is  in  fact  a  Cartesian  mapping  f sanction.  )   The 

imaginary  part  is  then  computed  as  the  conjugate  function, 
equation  (7)» 

In  addition  to  the  restrictions  on  the  velocity  dis- 
tribution mentioned  initially,  further  conditions  must 
be  met  in  this  method,  if  closed  contours  are  to  be  ob- 
tained.  Thus,  the  condition  for  closure  of  contour. 


n 


1/  '^ 


d^Z2  -  z^ 


-  1 


(M 
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and  the  recuired  coincidence  of   v„^   and   v^,   at 

infinity,  lead  to  the  follovi.'ing  three  restrictions  on 
the  real  part   R(cp)   of  the  integrand  In  equation  (3)+) 
considered  as  a  fionction  oi'   cp  in  the  circle  plane, 

JT2rr  '  /02tt  p,2v 

R(c)dcp=   (        R(cp)    COS    CD  dcp  =  R(C|))    sin   cp  dc,o  =  0         (55) 

0  JO  Jo 


Method   of   Wei nig   and    3ebeleln 

The  method  of  Weinig  and  Gebeleln  (reference  6)  nay- 
be  described  essentially  as  follovv's:  The  given  data  are 
the    same   as    in   the    3etz   method.      Consider   the   function 


V, 


log 


V 


^   =  log 

^1 


V. 


Vr 


-K^-  ^^1) 


{56) 


where   P„   and   P„    are  the  slopes  at  corresponding 
points  of  the  tv;o  airl'oils  (fig.  11).   Since  jvg^j   and 
Iv^ij   are  knov/n  functions  of   cp  with  the  data  as  given, 

and  since   log  — —      is  regular  outside  the  circle. 


V 


z- 


S^  _  -  f5^^   can  be  calculated  as  the  function  conjugate  to 

V, 


Z2    ^zi 


log 


V. 


a 


The  ane-le   P^    being  known,   S^_   is  thereby 

^        Z-i  o  J      ^2  "^ 


determined  and  hence,  by  simple  Integration,  the  unknov/n 
airfoil  coordinates  are  obtained. 

As  in  the  Betz  method,  the  condition  for  closure  of 
the  desired  contour 


'dz  =  '''  ^^'''^'^ 


ap 


Vr 


C  ^z 


dp  =  0 


(57) 


leads  to  the  additional  restrictions  on  the  prescribed 
velocitv  distribution. 
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_1_ 
2tt 


A 


2Tr 


log  Vg  ( c)|  dcp  =  0 


u 


0 


ri2Tr 


lOK 


V    (cp)|  sin  cp  dcp  =   -sin  2y 


> 


1     P"^^ 

-    I  log  V2(cp)|  COS    c?  dq^  =   -Tr(l    -   cos   2y) 

vJO 
v/here      y      is    given  by   equation    (5'3)« 


(38) 


Discussion  of  the  Various  Inverse  I'lethods 

The  methods  of  3etz  and  of  Weinig-Gebelein  may  be 
somewhat  narrower  in  scope  than  the  CMP  nethods.   The  use 
of  mapping  functions  such  as  in  equations  (53)  a-^d  (5^) 
Is  based  on  the  ability  to  specify   dZ2/dZ]_   unambiguously 

in  the  corresponding  regions.   This  requirement  appears 
to  restrict  the  contours  obtainable  by  these  methods  to 
those  bounding  simply  connected  regions.   Further  investi- 
gation of  this  point  is  necessary,  hov/ever.   By  the  CI'.IF 
methods,  figure-eight  contours  have  arisen  in  the  course 
of  solution  of  both  the  direct  and  the  inverse  problems. 
(See  the  g-per-cent  camber  derived  mean  line  (fig.  5(3-)  ) 
and  the  illustrative  examples  in  reference  1.)   Such  con- 
tours were  first  encountered  as  preliminary  results 
(unpublished)  in  using  the  method  of  potentials  v/ith  the 
Theodorsen-G-arrick  transformation.   The  CMF  apparently 
makes  no  fundamental  mathematical  distinction  between  simply 
connected  and  figure-eight  contours,  for  although  z    -    t, 
must  be  a  single -valued  function  of  z,     I,      or  p,   the 
coordinate   z   itself  is  of  the  same  char'acter  as  t      and 
the  latter  has  two  Riemann  slieets  at  its  disposal  in 
consequence  of  the  Joukowski  transformation  from  the  i- 
to  the  p-plane . 

The  methods  of  3etz  and  of  V/einig-Gebelein  require 
the  numerically  difficult  closure  conditions  fequations  (35) 
and  (58))  to  be  satisfied  in  adx^ance .   If  the  methods  are 
v/orked  through  for  prescribed  velocity  distributions 
v;hich  do  not  satisfy  these  conditions,  it  appears  that 
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open   contours   result.       In   the    ClIiF  methods,    however,    there 
is    either  no    closure    condition    (method   of   potentials)    or 
a   nuriierically   simple    Oxae    (method   of   derivatives): 

dcp 


J: 


dcp 


dC^  = 


3  ^^'  00 

[This  simple  closure  condition  in  the  method  of  deriva- 
tives is  fundamentally  a  consequence  of  the  fact  that 
the  required  absence  of  the  constant  term,  in  the  inverse 
pov/er  series  for  the  ClIP  derivative^  mapping  function 

(  ip  — — 7^ -,      mentioned  previously j  automatically  ex- 
cludes the  Inverse  first  pov/er  (the  residue  term)  from 
the  power  series  for   d(z  -  ^/dpfl   Thus,  physically 
impossible  velocity  distributions  lead  to  open  contours 
in  the  Betz-Welnig-Gebelein  methods  and  to  figure-eight 
contours  in  the  CMF  methods  (if  the  latter  converge). 
Prom  the  practical  point  of  viev;  in  these  cases,  it  may 
be  easier  to  obtain  the  airfoil  corresponding  to  the 
"best  possible"  -oh^^sicallv  attainable  velocitv  distri- 
but  ion  by  the  CI,IP  methods  than  by  the  others.   If  the 
succession  of  airfoils  deterxained  by  an  inverse  GI-IF  method 
is  seen  to  tend  toward  the  development  of  a  figure -eight, 
the  successive  approximations  can  be  stopped  at  the  "'best 
possible"  physically  real  airfoil. 

As  to  the  existence  and  uniqueness  of  a  solution  to 
the  inverse  problem  as  stated,  a  rigorous  discussion  of 
the  solutions,  for  a  prescribed  velocity  distribution, 
of  the  controlling  equations  (17)  ^  il'^)  j    snd  (9)  is 
lacking.   For  physically  possible  velocity  distributions, 
hovifever,  specified  as  a  function  of  percent  arc,  the 
Vifeinlg-Gebelein  method  shows  that  there  is  one  and  only 
one  airfoil  as  a  solution.   If,  hov;^ever,  the  velocity  is 
specified  as  a  function  of  percent  chord,  some  further 
condition  is  necessary.   This  requirement  is  evident  from 
the  fact  that  one  velocity  distribution  for  an  airfoil 
can,  for  differently  chosen  chords,  be  expressed  as  a 
different  function  of  percent  chord  in  each  case.   One 
chord  with  a  given  velocity  as  a  function  of  percent 
chord  can  therefore  have  more  than  one  corresponding 
airfoil.   There  is  reas-on  to  suppose  that  the  further 
condition  for  uniqueness  of  solution  in  this  case  is, 
the  chord  being  defined  as  in  the  section  "The  Direct 
Potential  Problem  for  Airfoils,"  that  the  ordinates  to 
the  airfoil  at  the  chordv^^ise  extrem^ities  be  specified. 
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From  the  experionce  with  the  CMP  methods  gained  to 
date,  it  is  believed  that  to  a  velocity  distribution 
specified  as  at  the  beginning  of  part  II,  and  v/lth  the 
further  condition  mentioned  in  the  percent-chord  case, 
there  corresponds  one  and  only  one  closed  contour  satis- 
fying the  CIvZP  system  of  equations.   It  is  furthermore 
believed  that  the  ClfF   nethods  are  flexible  enough  to 
converge  to  this  solution  in  at  least  those  cases  of 
aerodynamic  interest. 


CONCLUSIONS 


1.  The  conformal  transformation  of  an  airfoil  to 
a  straight  line  by   the  Cartesian  mapping  function  (CIIP) 
method  results  in  simpler  numerical  solutions  of  the 
direct  and  inverse  potential  problems  for  airfoils  than 
have  been  hitherto  available. 

2.  The  use  of  superposition  with  the  C?ff  method 
for  airfoils  provides  a  rigorous  counterpart  of  the 
approximate  methods  of  thin-airfoil  theory. 


Langley  Memorial  Aeronautical  Laboratory 

Rational  Advisory  ComLiittec  for  Aeronautics 
Langley  Field,  Va. 
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APPENDIX  A 

THE  CALCULATION  OP  CONJUGATE  FUNCTIONS 
EY  THE  RUWGE  SCHEDULE 


The  basic  calculation  for  the  type  of  mapping  func- 
tion treated  in  this  paper  and  in  reference  2  consists 
of  the  computation  of  the  real  part  of  an  analytic  func- 
tion on  a  circle,  given  the  imaginary  part,  or  vice 
versa.   To  this  end  the  conjugate  Fourier  series,  eojaa- 
tions  (i|.)  and  (5)j  or   the  conjugate  integral  relations, 
equations  (6)  and  (7),  ^^Q    available.   This  type  of  cal- 
culation appears  to  be  fundamental  in  many  kinds  of 
tv;o-dimenslonal  potential  problems.   For  example,  the 
solution  of  the  integral  equation  relating  normal  induced 
velocity  to  circulation  in  lifting-line  theory  can  be 
solved  easily  by  a  method  of  successive  approximation 
if  tho  transformation  from  the  ''lifting  line"  to  the 
circ.;.i-j  is  ^in.O't-Jn.      Quicker  methods  of  calculating  a  func- 
tion from  its  conjugate  than  those  given  in  this  appendix 
or  in  reference  2  v/ould  therefore  be  highly  useful. 

The  use  of  the  Fourier  series  rather  than  the 
integral  roj-ations  in  the  calculations  of  this  paper  was 
based  on  the  follov/ing  consideration.   Because  the  func- 
tion  l/z   is  regular  outside  the  unit  circle,  the  real 
and  imaginary  parts  of   l/z   on  the  unit  circle,  nam.ely, 
cos  cp   and   -sin  Q^,   satisfy  the  Integral  relations  (6), 
(7).   The  substitution  of   -sin  cp  for   A^r   in  equa- 
tion (6)  and  subsequent  numerical  evaluation  by  the  20- 
point  method  of   reference  2  gave  results  that  were  higher 
than  cos  <V     by  a  constant  error  of  2,8  percent.    Evalua- 
tion by  a  [|.0-point  method  reduced-  the'" error  by  half,  or 
to  l.[|.  percent.   By  the  Fourier  series,  on  the  other  hand, 
the  first  harmonic  (a  one-point  method)  suffices  to  give 
exact  results  in  this  case.   It  appears,  therefore,  that 
when  the  given  real  fianction  is  expressible  in  terms  of 
a  small  number  of  harmonics,  as  is  the  case  in  airfoil 
applications,  the  Fourier  series  method  is  preferable  to 
the  use  of  the  integral  relations. 

The  Runge  .schedule  offers  a  convenient  means  of 
carrying  out  the  basic  calculation  of  mapping  functions, 
namely,  the  analysis  of  a  periodic  function  into  its 
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Fouriei'  series  and  the  synohesis  of  a  Fourier  series 
into  a  function.   The  theory  and  use  of  the  schedule  is 
described,  for  example,  in  reference  S,  wherein  are  also 
given  schedules  for  12-,  2l\.~ ,    56-,  and  72-point  harmonic 
analyses . 

The  necessary  analyses  and  syntheses  in  the  direct 
and  inverse  CJ.IP  methods  are  carried  out  in  accordance 
with  equations  (li.)  and  (5)  9-nd  their  derivatives. 
Table  IX  contains  the  scheme  of  substitution  into  the 
Runge  schedule,  table  X,  for  the  various  CI'.CP  methods. 
In  the  direct  method,  for  example,  the  set  of  values 
5y/l2   corresponding  to  the  evenly  spaced  cp-values  is 
substituted  into  the   y^^   spaces  at  the  beginning  of  the 
sura-table.   The  sums  and  differences  of  these  quantities 
are  then  obtained  as  directed  at  the  left  of  the  indi- 
vidual tables  and  substituted  into  the  succeeding  tables. 
In  this  v/ay  the  entire  sum-table  is  filled  out.   Before 
the  product-table  is  used,  the  sum-table  should  be  checked. 

The  quantities  surrounded  by  the  heavy  lines  in  the 
sum-table  are  ne::t  multiplied  by  the  proper  factors  at 
the  left  of  the  product- table  and  the  results  entered 
in  the  appropriate  spaces  as  indicat^id  by  the  letters 
at  the  left  of  the  individual  product-spaces.   A  heavy 
horizontal  line  at  the  lov/er  left  edge  of  a  product- 
space  indicates  that  the  corresponding  prodiict  has 
already  been  obtained  in  a  previous  spsce  in  the  same 
rov/.   A  heavy  vertical  line  a?  ong  the  left  edge  of  a 
product-space  is  used  to  emphasize  that  tiie  negative 
value  of  the  product  of  the  sum-table  quantity  and  the 
product-table  factor  is  to  be  entered.   The  suras  of  the 
product-table  columns  are  then  entered  in  the  I,  II,  III, 
and  IV  spaces.   A  check  on  the  v/ork  of  the  product-table 
up  to  this  point  is  provided  by  the  columns  at  the  right. 
The  sujns  and  differences  of  the  I,  II,  III,  and  IV  quan- 
tities complete  the  product -table  and  give  the  Fourier 
coefficients   a,^,  bj^   corresponding  to   5y. 


In  order  to  perform  a  synthesis  calculation  from  a 
set  of  Fourier  coefficients   a^^,  h-^      to  the  values  of  the 
corresponding  function  at  the  even  cp-points,  the  coef- 
ficients  a^^,  bn   are  entered  in  the   d   and   D   spaces, 
respectively,  in  the  sum- table,  and  the  remainder  of  the 
simi-table  and  the  product-table  v/orked  through  as  before. 
The  final  values  in  the   a^-^,  bv^  spaces  of  the  product- 
table  are  then  entered  in  tl:e   d  and  D   spaces  at  thj.e 
beginning  of  the  sum-table  and  the  suras  and  differences 
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obtained  as  ind5.Gated  by  the  synthesis  column  at  the 
left.   (llote  that   dg   and  d^^^      ^^^    'to  be  multiplied 
by  2.)   The  resulting  y,-^   quantities  are  the  desired 
values  of  the  function. 

The  numerical  values  in  tables  X(a)  and  (b)  illus- 
trate the  process  of  obtaining  5x]_(cp)  from  &j\i^9)  in 
the  first  approximation  by  the  direct  CMF  method  for  the 
NACA  6512  airfoil. 
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APPENDIX  B 

Tl-ffi  MAPPING  OF  MORE  GEN-?3RAL  REGIONS 
Simply  Cormected  Regions 


If  the  CMP  nethod  Is  applied  to  the  napping  of  a 
simply  connected  boundary  v/ith  a  vertical  discontinuity, 
such  as  a  rectangle  or  an  infinite  line  with  a  vertical 
step,  the  ambiguity  of  the  ordinate   Ay   at  the  discon- 
tinuity will  prevent  an  automatic  and  rapid  convergence 
of  the  method.   Although  the  difficulty  could  be  lessened 
in  particular  cases  such  as  for  rectangles  by  taking  the 
diagonal  as  x-axis,  thus  removing  the  vertical  discon- 
tinuity, or  by  using  syr.iiTietry,  as  with  squares,  it  is 
evident  that  in  general  a  reference  shape  particularly 
suited  to  the  contour  under  investigation  is  needed. 
The  circle  has  been  shown  in  reference  2  to  be  a  good 
reference  shape  for  the  square.   It  could  be  expected 
therefore  that  an  ellipse  would  be  a  good  reference  shape 
for  the  rectangle.   i'Virthermore,  just  as  the  mapping 
function  based  on  the  circle  was  formed  of  an  anp^ular 
displacement  and  a  radial  displacement,  the  mapping 
function  based  on  the  ellipse  should  be  formed  of  dis- 
placements along  and  orthogonal  to  the  ellipse,  that  is, 
should  be  specified  by  elliptic  coordinates.   The  speci- 
fication of  a  figure  by  elliptic   coordinates  (^ ,    6)   in 
the  physical  plane   z   is  equivalent,  however,  to  the 
transformation  of  the  figure  to  a  t' -plane  by  the  two 
transformations 


z  =  p'  + 


P 


where   p' 


\i/+i0 


t'  =  log  p'   where   t'  =   ^\t  +    19 


(59) 


and  specifying  the  transformed  figure  by  the  Cartesian 
coordinates  of  the  t ' -plane  {^l^ ,    9).   The  rectangle  under 
consideration  v/ill  be  a  near-circular  shape  in  the  p*- 
plane  and  a  near-straight  line  shape  in  the  t' -plane. 
The  mapping  of  the  rectangle  by  means  of  an  elliptic 
mapping  function  in  the  physical  plane  is  therefore  seen 
to  be  accomplished  by  the  Theodorsen-Garrick  method  in 
the  near-circle  p' -plane  and  by  the  CMF  method  in  the 
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near-stralghb  line  t' -plane.   Prom  this  point  of  vlevi/, 
therefore,  the  Theodorsen-Garriek  method  consists  of 
specifying  an  airfoil  in  the  physical  plane  by  elliptic 
coordinates,  forming  the  corresponding  elliptic  mapping 
function    ( i|/  -  xif^ )  -  i  e ,   which  conformally  relates  the 
airfoil  to  an  ellipse  or  Joukovvski  airfoil  as  a  basic 
shape,  and  expressing'  the  elliptic  mapping  function  as  a 
regular  function  outside  the  circle.   On  the  other  hand, 
in  the   t'  =  log  p' -plane  the  Theodorsen-G-arrick  method 
consists  of  the  transformation  of  the  near-straight 
line  ii{Q)      to  the  straight  line   ^^  =  Constant  by  means 
of  what  is  now  the  GMP  -  U'   ~   ^q)    -    ie.   Thus,  the 
Theodorsen-Garrick  method  may  be  regarded  as  a  form  of 
the  CRIP  method,  in  which   log  p'   takes  the  place  of   z 
and  log  p,   the  place  of   t. 

The  mapping  of  simply  connected  regions  by  dif- 
ference mapping  functions  based  on  the  curvilinear  co- 
ordinates appi^opriate  to  the  particvilar  reference  shape 
considered  is  therefore  equivalent  to  using  the  GMF  dif- 
ference function  z    ~    I      in  the  plane  of  the  near-straight 
line  into  which  the  reference  shape  is  initially  trans- 
formed. 


Mapping  of  the  Entire  Field 

The  Fourier  series  representation  of  mapping  func- 
tions, equations  (ii)  and  (5)*  enables  the  calculation  of 
corresponding  points  in  the  tvs/o  regions  to  be  made,  once 
the  correspondence  of  the  boundaries  has  been  calculated. 
By  the  latter  calculation  the  coefficients   aj^,  h-^-^      and 
the  radius   R   of  the  circle  of  correspondence  have  been 
determined.   If  now  a  large:.^  radius   R'  >  R   be  substi- 
tuted for   R   in  equations  ([].)  and  (5)^  the  resulting 
synthesis  of  the  Fourier  series  will  yield  the  mapping 
function  for  the  circle  of  radius   R' ;   that  is,  v/ill 
determine  points  in  the  given  plane  corresponding  to  the 
points  in  the  circle  plane  at  the  distance   R'   from  the 
origin.   It  is  nocessary,  of  course,  to  use  the  mapping 
function  in  conjunction  v^ith  the  shape  in  the  physical 
plane  corresponding  to  the  larger  circle.   In  this  way  the 
entire  corresponding  fields  can  be  mapped  out.   It  m.ay 
be  noted  that  substitution  of   R'  <  R   for   R   in  equa- 
tions (L.)  and  (5)  enables  the  mapping  of  those  corre- 
sponding points  inside  the  original  contours  for  v/hich 
the  resulting  Fourier  series  converge. 
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It  appears  to  be  r.oro  C'ifficult  to  find  the  point 
in  the  circle  plane  corresponding  to  a  point  of  the  given 
plane  than  vice  versa.   This  calculation  may,  hov/ever, 
be  accomplished  by  a  nethod  of  successive  approximations. 
P'or  example,  if  the  givexi  plane  is  that  of  a  near  circle 
the  polar  coordinates  of  tlie  given  point  in  the  near- 
circle  plane  are  assumed  to  be  a  first  approximation  to 
the  coordinates   R'   and   cp  of  the  desired  point  in  the 
circle  plane.   Substitution  of  these  values  into  equa- 
tions (Ij.)  and  (5)  yields  a  first  approximate  mapping 
function  which  can  be  used  to  correct  the  coordinates   R' 
and   CO,  etc  . 


Biplanes 

In  the  case  of  the  biplane  arrangement  the  CLIP  may 
be  set  up  directly  in  the  physical  plane  in  the  same  v/ay 
as  for  the  single  airfoil.   In  place  of  the  simple  trans- 
formation from  straight  line  to  circle,  however,  the 
transformation  from  the  two   extended  chord  lines  of  the 
airfoils  to  two  concentr'lc  circles  Is  used.   This  trans- 
formation is  derived  in  reference  9*   The  CMP  method  for 
biplanes  bears  the  same  relation  to  the  method,  of  ref- 
erence 9  that  the  CMP  method  for  monoplane  airfoils  bears 
to  the  Theodorsen-Garrick  method  (reference  2). 

For  biplanes  (fig.  12)  the  CMF   z  -  t,  being  regular 
in  the  region  outside  the  two  straight  lines,  is  regular 
in  the  annular  region  of  the  p-plane  and  consequently  is 
expressible  as  a  Laurent  series  in  p 


^=> 


CO  \^ 

pn 


v>;here 


(1^0) 


°n  =  ^n  ^  i^n 

If,  for  the  inner  circle,  the  relationship  is  written 

z  -  ^  =  Ax-[_  +  i  Ay-^ 

(liD 
p  =   R-^e 


ko 
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and  for   the    outer   circle 


z    -   l    =  Ax^   +   i   A 


p   =   RpO 


icp 


72 


> 


(142) 


there  Is  obtaD.ned,  upon  substitution  into  equation  ( Lj-O ) 
and  reduction 


Ax-i  (q)  =  s.^  + 


^n^^-n 


R- 


n 


cos   ncp 


A 


^n  -  ^-n 


-n  ncp      (Ij-Ja) 


Ax,(cp)  =  a     +\        5a^-^^ 


R, 


n 


cos   np+\        — -^^  sin  ncp      ik-^h) 


1 


Ay-j_(cp)  =bo 


-^y?. 


cp 


=  bo  + 


1 


l^n+b-n 


■p. 


n 


cos   ncp  -  \ 


^n  ~  ^--n 


R 


n 


2 


sin  nCP      (i^Jd) 


These  equations  are  the  generalization  to  the  biplane  of 
equations  (I4)  s.nd  (5)'  I'be  corresponding  integral  rela- 
tions may  be  derived  as  In  reference  9 • 


The  solution  of 
or  the  inverse  prob] 
successive  approxiria 
method  the  tvro  alrfo 
mat ion  biplane  were 
be  taken  as  the  inlt 
formation  of  referen 
tions  on  the  straigh 
evenly  spaced   Q  po 


equations  (Ij-3)  Iri  either  the  direct 
em.  nay  be  accomplished  as  before  by 
tions.   For  exaraple,  in  the  dii'^ect 
lis  are  given.   If  no  initial  approxi' 
available,  the  tvjo  chord  lines  would 
lal  straight  lines.   By  the  trans- 
ce  9  this  fixes  the  chordv/ise  loca- 
t  lines  corresponding  to  a  set  of 
ints  on  the  concentric  circles.   The 
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ordlnates  i\ji^{fp)      can  therefore  be  ineasured,  v/hlch 
determines   Ayp(cn)   by  analysis  and  synthesis  ox"'  equa- 
tions (i|5o)  and  (lijd),  respectively.   (The  radius  ratio 
R2/R]L   is  fixed  by  the  initial  transformation  frori  the 
straight  lines  to  the  concentric  circles.)   These   Ay2((p) 
values  then  determine  a  set  of      LXpi'V)    values  by  the 

given  shape  of  the  second  airfoil  and  the  known  chordv/ise 
locations  of  its  first  appi-^oximation  straight  line. 
Analysis  of   Ax^  (<P )   and  subsequent  synthesis  of  l\x-,{'V) 

by  equations  (i4.5b)  and  (i^a),  respectively,  determines  a 
correction  to   R-,   by  a  horizontal  stretching  process 

(constant   Ax,  Ay  -  adjustr.ient  of   r-,  )  to  riaintain  the 
given  airfoil  chord.   The  procedure  is  now  repeated  v/lth 
^Jpi''^)      as  the  initial  set  of  laeasuri^id  ordinates  that 
determines   Ay^(Ci:'),  t:X^{<V) ,      and  Ax^(cp)   as  before.   The 

radius   R2   can  now  be  similarly  corrected.   This  step 
completes  the  first  approximation.   For  the  second  approxl- 
m.ation  a  new  corresnondence  betv/een  the  corrected  straight 
lines  and  the  concentric  circles  is  calculated  and  the 
procedure  repeated. 

The  Inverse  problem  could  also  be  solved  by  m.ethods 
similar  to  those  given  for  the  isolated  airfoil.   Sup- 
pose, for  examole,  a  vdng  section  were  given  and  it  were 
desired  to  derive  a  slat  of  given  chord  and  given  approxi- 
mxate  location  and  having  a  prescribed  velocity  distri- 
bution.  The  method  of  surface  potentials,  for  example, 
enables  the  calculation  of  a  first  approximate   Ax-]_(<p) 

(subscript  1  refers  to  slat).   The  initial  correspondence 
of  points  between  the  straight  lines  and  concentric 
circles,  and  therefore  also   R2/R1,   being  determined  by 
the  initially  assumed  straight  lines,  the  function  Ax2(cp) 
is  thereupon  obtained  by  analysis  and  synthesis  of  equa- 
tions (I|-3a)  and  (kj'o)*  respectively.   The  horizontal  dis- 
placem_ent   Ax2  (  9)   thence  determines   Ay2(C))   by  the 
known  shape  of  the  main  wing  section.   The  determination 
of   Ay2_(cD)   by  analysis  and  synthesis  of  equations  (Ij-Jd) 
and  (I|.5c)  completes  the  calculation  of  the  first  approxl- 
m^ate  slat  section,  for  whicli  the  exact  velocity  distri- 
bution can  now  also  be  calculated.   If  the  main  ?/ing 
section  were  also  unknov/n  then  the  Vi/ing  section  above  Is 
regarded  as  an  Initial  approximation,  the  role  of  the  two 
airfoils  is  reversed,  and  the  procedure  repeated  to  com.- 
plete  the  first  approximation. 
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The  Clip  •metlicd  can  bs  generalized  In  the  same  manner 
for  multiply  connected  regions.   The  transformation  from 
the  n  reference  shapes  (such  as  straight  lines)   to  n 
circles  being  presumed  known,  the  CMP  can  be  set  up  as  a 
series  convergent  in  the  region  between  the   n   circles, 
and  the  mapping  function  for  each  boundary  explicitly 
expressed  by  allov/ing  the  coordinate  vector  to  assume  its 
value  on  each  bou.ndary  in  turn.   A  method  of  successive 
approximation  for  the  solution  of  the  resulting  equations 
depending  on  the  particular  problem  under  consideration 
would  then  be  established. 


Cascade  of  Airfoils 

A  simplified  but  practically  Important  n-body  problem, 
namely,  the  cascade  of  airfoils,  may  be  mentioned  finally. 

The  reference  shape  into  which  the  cascade  of  air- 
foils, figure  15,  is  to  bo  transformed  is  chosen  as  the 
cascade  of  strlght  lines  coinciding  v/ith  the  extended 
chord  lines  of  the  airfoils  of  the  cascade.   The  trans- 
formation from  the  cascade  of  straight  lines  to  a  single 
circle  Is  well-known,  reference  10.   The  C¥F   chosen  as 
indicated  in  figure  I5  is  therefore  expressible  as  an 
inverse  power  series  In  tl^e  circle  plane  and  the  resulting 
procedure  in  either  the  direct  or  the  inverse  problem  is 
seen  to  be  essentially  the  same  as  for  Isolated  airfoils. 
The  detailed  application  of  the  CLIP  to  cascades  of  air- 
foils Is  given  in  reference  1. 
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APPEITDIX   C 

THE  DETERMINATION  0?  ivJ/.PPING  FUNCTIONS  BY  THE 
CAUCIiY  INTEGRAL  FORMULA 

The  foregoing  nethods  of  conformal  transformation 
heve  been  presented  from  the  point  of  viev;  of  represen- 
tation of  the  various  mapping  functions  as  infinite 
series.   In  particular,  the  expression  of  the  Cartesian 
mapping  function  as  an  inverse  power  series  valid  every- 
where outside  and  on   a  circle  led  to  the  Fourier  series 
representation  for  the  CMP  on  the  circle  itself.   The 
integral  formula  representation  vvas  then  obtained  from 
the  Fourier  series  hj   the  method  of  reference  5«   It  is 
of  interest  to  see  hov/  the  integral  i-elations  (6)  and  (7) 
can  be  derived  directly  from  the  Cauchy  integral  formula 
for  a  function  regular  outside  a  circle.   (These  integral 
relations  have  also  been  derived  by  Betz,  reference  7, 
by  a  hydrodynamical  argument.)   Since  the  applicability 
of  the  Cauchy  integral  formula  is  not  restricted  to 
circular  boundaries,  hovifever,  the  results  v/ill  be  capable 
of  generalization,  in  principle  at  least,  to  arbitrary 
simply  and  multiply  connected  regions. 

The  Cauchy  integral  formula  gives  the  values  of  an 
analytic  function  f  (p)   v/ithin  a  simply  connected  do- 
main D   in  terms  of  its  values   f(t)   on  the  boundary 
of  the  domain  as 

f(p)  =  _k_.  r£iLL  dt  ikh) 

dT(±   J   t  -  p 

v/here  the  path  of  integration  is  counterclockwise  around 
the  boundary.   Consider  the  domain   D   outside  the  sim.ple 
closed  boundary  C   in  the  p-plane  (fig.lJ4).   This  domain 
can  be  made  simply  connected  by  an  outer  boundary  B   and 
the  cuts  between  the  two  boundaries,  as  indicated  by  the 
dotted  lines.   The  Cauchy  integral  formula  for  the  func- 
tion f(p)   at  an  interior  point   p   of  the  domain  D, 
in  terms  of  its  values  on  the  boundary  is 


f  (p)  =  _L.  f  JUL  ,t  *  -J- 

2Tri     I      t    -   p  2rfi 


—   at  (45) 

t    -    p 
B 
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where  the  equal  and  opposite  Integrals  along  the  cuts 
have  been  oinittod.   "^'he  paths  of  integration  are  indicated 
b^""  the  arrows  in  figure  llj..   The  function  f (p)   is  as- 
sumed to  be  regular  everywhere  outside  the  boundary  C 
and  in  particular  to  approach  the  limiting  value  t^   as 
p — ^<».   If  the  boundary   B   is  enlarged  indefinitely, 
the  Integrand  of  the  second  integral  of  equation  (14-5) 
approaches   fco/^  ^^^   thus 


f(t) 
t  -  p 


dt  =  f^  (i|6) 


The  variable   p  v/ill  now  be  made  to  approach  a  point   t' 
on  the  boundary   C,   and  equation  (I4-5)  will  consequently 
reduce  to  an  integral  equation  for  the  boundary  values 
of  a  function  regular  ever^rv/here  outside  and  on  the  boun- 
dary.  In  order  to  evaluate  properly  the  contribution  of 
the  remaining  (first)  integral  of  equation  ( )-i-5 )  ^^   "ti^e 
neighborhood  of   t',   the  boundary  G   is  modified  as 
indicated  in  figure  ll^.   The  point  p   is  made  the  center 
of  a  semicircle  whose  ends  are  faired  into  the  original 
boiindary.   As   p — jt',   the  modified  boundary  approaches 
coincidence  with  the  original  boundary.   The  integral 
over  the  modified  boLindary  is  novv  evaluated  as  the  sum 
of  the  Integral  over  the  semicircle,  which  in  the  limit 

is  half  the  residue  of  the  integrand  or   —  f(t')^  and  the 

integral  over  the  rest  of  the  path,  v/hich  in  the  limit 
is  analogous  to  the  Cauchy  principal  value  of  a  real 
definite  integral  of  vuhich  the  integrand  becomes  infinite 
et  som.e  point  in  the  interval  of  integration.   Equa- 
tion ( I1.5 )  therefore  becomes,  in  the  limit. 

In  addition,  there  is  the  auxiliary  condition  that 

(ii8) 


2tt1  I 
u 


vjhich  folloviis  from  the  fact  that   f  (p)   is  regular  every- 
where outside  the  boundary   0.   Equation  (IlJ)    is  well 
knovm  in  the  theory  of  functions  of  a  complex  variable. 
(See  reference  11.) 


If,  no"^^ ,    th'3  function  f(p)   is  taken  as  the  Cartesian 
napping  function    z    -    l        or,  on  the  boundary, 

f(t)  =  Ax  +  i  Ay  ();9) 

and  if,  further,  the  bounda?'';^,'   C   is  taken  as  a  circle 
v/ith  origin  at  the  center. 


/■o+W 


1- .'  =  A  ^ 


r  (5-) 


substitution  of  equations  ([i9)  and  (50)  into  equation  (lj-7) 
and  using  equation  (I4.8)  (v/ith  f^  =  0)  leads  to  the  inte- 
gral relations  (6)  and  (7)-   I-t"  the  polar  mapping  func- 

tion   log  =:^  =  'I'  -  ie  E  (\j/  -  \J/^)  -  l(cp  -  9)   (reference  2) 

is  substitiited  for   f(t),   the  Theodorsen-Garrick  integral 
relations  are  obtained. 

The  Cauch^-  integral  formula  has  already  been  applied 
(reference  12)  to  problems  of  conformal  mapping  in  the 
manner  just  indicated.   Bergman  has  included  in  refer- 
ence 12  (chapter  ZI )  contributions  of  two  Russian  authors, 
Gershgorln  and  Krylov.   In  reference  12  the  mapping  func- 
tion from  a  circle  to  a  near  circle  viras  taken  in  the  form 
log  p.   The  resulting  integi'al  equation  does  not  appear 
to  be  as  convenient  as  those  of  the  C^iP  methods .   The  use 

of  forms  such  as   log  - —  or  z    -•    i      are  not  onlv  accurate 

p 

nuTiierically  since  they  express  changes  in  the  coordinates 

of  the  boundaries,  but  also  they  lead  to  pairs  of  integral 

eo^uations  v/hlch  contain  the.  solutions  of  both  the  direct 

and  the  inverse  potential  problems . 

Prom  the  analysis  given  It  appears  possible  to  trans- 
form conforraally  from  one  boundary  to  another  vifithout 
requiring  the  transformation  from  either  boundary  to  a 
circle,  since  the  boundary   G   In  equation  (l|-7)  can  be 
rather  arbitrary  and  f (t)   can  be  taken  as  a  mapping 
function  from  this  boundary  to  another  arbitrary  one. 
The  resulting  integral  equation  for  the  mapping  function 
Is,  hov/ever,  not  as  easy  to  solve  numer^ically  as  when  the 
boundary  C   is  a  circle. 
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Once  the  confor'nal  corresponde-nce  between  two  boun- 
daries Is  knovv'n,  corresponding  points  outside  the  boun- 
daries can  be  determined  by  the  Cauchy  integral  formula 
(Lii).).   It  is  noted  that  the  Cauchy  Integral  gives  the 
correspondence  of  Individual  pairs  of  points  rather  than 
the  correspondence  of  entire  boundaries  at  once,  v/hich 
is  given  by  the  Fourier  series  representation.   Further- 
more, the  possibility  exists  of  determining  pairs  of 
corresponding  points  Inside  the  given  boundaries  by  the 
Cauchy  integral,  that  is,  of  analytically  continuing  the 
conformal  transformation  beyond  the  original  domains. 
For  if  the  transformation  from  a  boxindary  C   in  a 
p-plane  to  a  boiindary  C'   in  a  p' -plane  vi?ere  knovm,  the 
outside  regions  corresponding,  then  the  correspondence 
between  a  boundary  A   internal  to   C   and  a  boundary  A' 
internal  to   C',   if  It  existed,  could  be  determined  by 
an  application  of  the  Cauchy  integral  formula  to  the 
region  bounded  by  A   and  C. 

For  example,  if  the  boundaries   A   and   C   are  taken 
as  concentric  circles  and  the  mapping  function  as 

f (p)  =   log  ^ 

=  ^l/  -  1^  =  (il/  -  X)  -  Kcp  -  0)  (51) 

in  the  notation  of  figure  15,  the  Cauchy  Integral  formula 
applied  to  the  annular  region  in  the  p-plane  (assum.ed 
free  of  singularities  of  the  mapping  function)  yields, 
in  the  limit  as  the  variable  point   p   approaches  the 
inner  circle  A, 

2"Jo  cosh  (Xq  - Xi)  -  cos  (^^  -  CP^') 

1    p2TT    .  Cpi   -   C?    • 

^l(^^l')=-   27  '^^^l)    cot  -i-^:-J- dcpi 


d'Po  (52a) 


2Mo  cosh  (X,-Xi)- cos   (P,  -CPi') 
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In  adoitlon,  the  condition  of  regularity  of  the  function 
f(p)  in  the  annular  region  yields  the  auxiliary  condi- 
tions 


1    r2^        ^  1      />2TT 

2tt  /  1      1        2rr    /  ^    ^' 


Jo 


1      ,^ 


1      P2TT 


2Tr 


^o^^o 


55) 


In  the  problem  under  consideration,  the  mapping  function 


*o<^'o'  - 


1S<%' 


for  the  cuter  boundaries 
of 


s  lea  own. 


-o 


2^1 


the  concentric  circles   are  given, 
of  equations  (52)  are  thus  known  functions  of  •?]_' 


The  radii   e'^^,   e 

The  second  integrals 
jc    c^xw  ^...,^    ..x..^.,..  ^..ux.^^^.^.x^  ^^       .  j_  .   Equa- 
tions (52a)  and  (52b)  therefore  constitute  a  pair  of 
integral  equations,  similar  to  those  of  Theodorsen-Garrick, 
for  the  mapping  function  '?^i(''?i)  -  ^^if'^iV   pertaining 
to  the  Irxner  boundaries.     ^   ■       ^   •' 


It  is  noted  that  if  the  variable  point   p   of  the 
Cauchy  integral  formula  for  the  annular  region  is  made 
to  a^'proach  the  outer  boundary  G,   then  two  additional 
integral  equations  similar  to  equations  (52a)  and  (52b) 
are  obtained.   These  equations,  together  with  equa- 
tions (55),  are  a  generalization  to  the  case  of  ring 
of  the  corresponding  Theodorsen-Garrick 


equations  for  slmpljr  connected  regions  and 
for  the  conformal  mapping  of  near  circular 


can  be  used 
ring  regions. 
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30- PERCENT   THICI<i\ESS    JOUKOV/SKI    PROFILE 
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(radians ) 



^Jo 

dAxo/dcp 

dAy^/dcp 

^      0    X     ^ 

12 

1 

-0.519 

-.301^ 

0 
.096k 

0 
.119 

0.575 

.555 

2 

-.258 

.132 

,226 

.296 

5 

-.190 

,250 

.309 

.206 

k 

-.101 

.287 

.352 

.089I1- 

5 

- .  0072i| 

.295 

.352 

-.0551 

6 

.0798 

.270 

.301^ 

-.11^9 

7 

.11^8 

.213 

.212 

-.233 

8 

.189 

.150 

.0916 

-.272 

9 

.197 

.0810 

-.0261 

-.21+0 

10 

.179 

.0291 

-.0958 

-.11+9 

" 

.15^ 

.00[!.12 

-.082k 

-.031+6 

1 12 

.i)i? 

0 

0 

0 

TABLE    II 
CONST/iNTS   USED  WITH      CMF      OF   TABLE    I 


Profile 

T 

Xu^ 

T 

r 

cPt 
(d.eg) 

^t 

Joukowski 

0.30 

1.000 

0.0887 

1.230 

0 

180 

1.000 

Derived 

.21+ 

.805 

.  0716 

1.185 

0 

180 

.835 

Derived 

.12 

.1+02 

.0357 

1.0928  0 

180 

.I1-53 
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TABL3    III 
CMF      FOR   6-PEJ?CENT-CAMBER   CIRCULAR-ARC    PROFILE 


(radians) 

Ax, 

^^0 

dAx^/dcp' 

dAy^/d?) 

6    X  -^- 

12 

0 

0,120 

0.108 

0 

7 

.0270 

.111+ 

.0960 

- .  0^81| 

8 

.014.82 

.0955 

.0638 

-.0871 

9 

.0592 

.0691+ 

.0171 

-.109 

10 

.0565 

.OJ4O5 

-.0563 

-.106 

11 

.01^08 

.0160 

-.O8I1I+ 

-.0781 

12 

.011x2 

.00169 

-.115 

-.0279 

13 

-.0170 

.002l|6 

-.117 

.031^6 

ik- 

-.ol|59 

.0191I 

-.0852 

.0926 

15 

-.0587 

.014.90 

-.0239 

.128 

i6 

-.0552 

.0828 

.0506 

.123 

17 

-.0535 

.110 

.113 

.0756 

18 

0 

.120 

.136 

0 

TABLE    IV 
CONSTANTS   USED  VVITH      CMF      OF      TABLE    III 


Profile 

r 

(per- 
cent ) 

r—— 

^c 

Xu^ 

T 

r 

(deg) 

Cpr[n 

(deg) 

(deg) 

^Udeal 

Derived 

-7 

9 

0.502 

0.501 

0 

1.0052 

-5.37 

183.37 

0 

0.37 

Circular 
arc 

6 

1.000 

1.000 

0 

1.0072 

-6.81^ 

186.81^ 

0 

.75 

Derived 

9 

1.502 

1.1^99 

0 

1.0050 

-10.26 

190.26 

0 
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TABLE   IX 
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